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I. INTRODUCTION 


If modern aircraft antennas are to function properly, the antenna 
patterns must meet certain system requirements. The conventional design 
procedure for an antenna on a particular aircraft has been to evaluate 
the performance of a candidate antenna system based on numerous 
scale-model measurements. This approach not only requires a great deal 
of engineering time and expense, it also has the following drawback. 

When the near-field patterns are relatively easy to measure, the 
far-field patterns are not. To obtain the far-field patterns 
cost-effectively, there has been a great deal of interest in determining 
far-field patterns based on near-field pattern measurements. Currently, 
most of the attention has focused on plane, cylindrical and spherical 
wave spectrum approaches. However, each of these spectrum approaches is 
basically an integral transform which in itself can be tedious and 
expensive. Thus, the following dilemma prevails: far-field patterns 

are desired but cannot be easily measured directly; near-field patterns 
are much easier to measure but cannot be simply transformed to the far 
field. 


One approach to solve the above problem is the geometrical theory 
of diffraction (GTD). It is a high frequency technique which is valid 
for computing both the near-field and far-field patterns when the source 
and various scattering centers are separated on the order of a 
wavelength. The near-field solution can be easily verified by a 
near-field measurement. Once this near-field verification is 
accomplished, the solution can be directly extended to the far field 
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without the need of a transformation. The GTD is directly applicable to 
the near-field pattern prediction because the GTD postulates are not 
violated in the sense that the receiver is essentially in the far field 
of each isolated specular point. For instance, the receiver might not 
be in the far field of a flat plate; yet, it is sufficiently removed 
from each of the edge diffraction points that the GTD is valid (i.e., 
the receiver is at least a wavelength away from the isolated diffraction 
points). Consequently, a GTD solution can effectively solve both the 
near-field and the far-field airborne antenna patterns. 

This Fortran 77 computer code has been developed at Ohio State 
University to investigate the radiation patterns of antennas mounted on 
an aircraft fuselage which is modeled by a composite ellipsoid. The 
computer code is used to compute the near zone and far zone radiated 
fields for antennas mounted on a composite ellipsoid and in the presence 
of a set of finite flat plates. The analysis applied in the development 
of this code is based on the uniform geometrical theory of diffraction 
(DTD) [1,2,3]. 

The code allows the user to simulate a wide variety of complex 
electromagnetic radiation problems using the ellipsoid/plates model. 

For example, the composite ellipsoid can be used to accurately simulate 
the fuselage of an aircraft; whereas, the plates are used to represent 
the wings, stabilizers, stores, etc. Alternatively, the antenna could 
be mounted directly on a ship mast. In this case the mast could be 
approximated by the composite ellipsoid with the other ship structures 
simulated by flat plates. Note that the plates can be attached to the 
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composite ellipsoid and/or to other plates. In fact, the plates can be 
connected together to form a box. This code is specifically designed to 
analyze the radiation characteristics of antennas mounted on aircraft 
configurations. 

As with any ray optical solution such as this UTD code, there is a 
limit to the number of interactions included in the field computation. 

In this case, the code includes the source, reflected, diffracted, and 
higher order terms such as the reflected/reflected , 
reflected/diffracted, diffracted/reflected, and diffracted/diffracted 
fields. The higher order terms are due to the multiple field 
interactions between the simulation plates. It assumes that the 
higher-order diffracted and reflected fields from the composite 
ellipsoid surface are small and can be neglected. The user may request 
the code (by using the “TO:" COMMAND) to compute the higher order terms 
when he thinks they have a significant effect on the results; otherwise, 
the code will compute first order terms only. This implies that the 
code can handle structures for which the energy does not significantly 
bounce back-and-forth across the target. In any event, the code 
automatically shadows all terms, such that if a higher-order interaction 
should have been included the resulting pattern will contain a 
discontinuity. These higher-order terms are normally negligible and can 
only affect the pattern in rather small sectors. However, if they are 
significant in some region, the amplitude of the jump is associated with 
the radiation level of the missing higher-order term. Consequently when 
the solution fails because of a lack of higher-order terms, it tends 
to indicate its failure. 
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The code has the flexibility to handle arbitrary pattern cuts. In 
addition, an arbitrary antenna type can be analyzed provided the current 
distribution across the aperture is known. This is done by 
approximating the distribution by a set of magnetic current elements 
mounted on or electric currents normal to the composite ellipsoid 
surface. The magnetic current elements have a cosine distribution 
along the magnetic current direction and a uniform distribution in the 
orthogonal direction. The normal electric current represents a monopole 
provided that its length isn't greater than a quarter wavelength. 

The mutual coupling effect for monopole arrays mounted on a 
fuselage can be handled by thin-wire theory [41, if the region near the 
array is nearly flat. For engineering purposes, image theory can be 
applied to calculate the relative current distributions as equivalent 
dipole arrays. The relative current value on each dipole is then taken 
to be part of the input data for each monopole source specification. 

The final pattern is the superposition of the contributions from each 
individual monopole. 

The limitations associated with the computer code result from the 
basic nature of the analyses. The solution is derived using the UTD, 
which is a high frequency approach. In terms of the scattering from 
plate structures this means that each plate should have edges at least a 
wavelength long. In terms of the composite ellipsoid structure its 
major and minor radii should be at least a wavelength in extent. In 
addition, each antenna element should be at least a wavelength from all 
edges. In some cases, the wavelength limit can be reduced to a quarter 
wavelength for engineering purposes. 
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The present code requires approximately 707K bytes of storage. It 
will run a pattern cut of 360 points for a commercial aircraft model 
(Example 3, 6 plates included) with one antenna element in approximately 
4 minutes on a VAX 11/780 Computer. 

This user's manual is designed to give an overall view of the 
operation of the computer code, to instruct a user in how to use it to 
model structures, and to show the validity of the code by comparing 
various computed results against measured data whenever available. 
Section II describes an overall view of the organization of the program. 
The definition of the input is given in Section III. How to apply the 
capabilities of this input data to a practical structure is briefly 
discussed in Section IV. This includes a clarification of the subtle 
points of interpreting the input data. The representation of the output 
is discussed in Section V. Numerous practical airborne antenna problems 
are presented in Section VI to illustrate the operation, versatility, 
and validity of the code. 

II. PRINCIPLES OF OPERATION 

The analytical modeling of complex scattering shapes in order to 
predict the radiation patterns of antennas has been accomplished using 
the Uniform Geometrical Theory of Diffraction (UTD) [1,2,3], This is a 
high frequency technique that allows a complicated structure to be 
approximated by basic shapes representing canonical problems in the UTD. 
These shapes include flat or curved wedges and convex curved surfaces. 
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The UTD is a ray optical technique, and it, therefore, allows one to 
gain some physical insight into the various scattering and diffraction 
mechanisms involved. Consequently, one is able to quickly seek out the 
dominant mechanisms for a given geometrical configuration and radiation 
sector. This, in turn, leads to an accurate engineering solution to 
practical antenna problems. This approach has been used successfully in 
the past to model aircraft shapes [5,6,7,8,9,10] and ship-like 
structures [11,12,13]. 

This section briefly describes the basic operation of this code for 
the analysis of antennas in an aircraft environment. The present 
version of the code allows the analysis of objects that can be modeled 
by flat plates and a composite ellipsoid all of which are built up from 
the basic canonical problems. These shapes allow one to model a wide 
variety of structures in the UHF range and above where the scattering 
structures are large in terms of a wavelength. The general rule is that 
the lower frequency limit of this solution is dictated by the spacings 
between the various scattering centers and their overall size. In 
practice this means that the smallest dimensions should be on the order 
of a wavelength. 

The positive time convention ejut has been used in this code, and, 
all the structures are assumed to be perfectly conducting and surrounded 
by free space. 

As mentioned above, the UTD approach is ideal for a general high 
frequency study of aircraft antennas in that only the most basic 
structural features of an otherwise very complicated structure need to 
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be modeled. This is because ray optical techniques are used to 
determine components of the field incident on and diffracted by various 
structures. Components of the diffracted fields are found using the UTD 
solutions in terms of the individual rays which are summed with the 
geometrical optics terms at the field point. The rays from a given 
scatterer tend to interact with other structures causing various 
higher-order terms. In this way one can trace out the various possible 
combinations of rays that interact between scatterers and include only 
the dominant terms in the solution. Thus, one need only be concerned 
with the important scattering components and neglect all other 
higher-order terms. This method leads to accurate and efficient 
computer codes that can be systematically written and tested. 

Complex problems are built up from similar components in terms of a 
modular computer code. This modular approach is illustrated in the 
block diagram of the main program shown in Table I. The code is broken 
up into many subroutines that represent different scattered field 
components, ray tracing sections, geodesic path algorithms [14,151, 
plate attachment calculation and shadowing routines. One is referred to 
Reference [161 for more details on this topic. As can be seen from the 
flow chart, the code is structured so that all of one type of scattered 
field is computed at one time for the complete pattern cut so that the 
amount of core swapping is minimized, thereby, reducing overlaying and 
increasing efficiency. The results are then, superimposed in the main 
program as the various segments are executed. 
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TABLE I 

BLOCK DIAGRAM OF THE MAIN PROGRAM 


r 

A 

RW 


V. 

J 




ECHO INPUT DATA 
SEE TABLE II) 


COMPLETE 




DEFINE PATTERN COORDINATE 
INFORMATION 

DETERMINE ALL FIXED GEOMETRY 
INITIALIZE TOTAL FIELD TO ZERO 


COMPUTE VARIOUS UTD TERMS 
(NOTE: PATTERN LOOP IN EACH UTD 

TERM SUBROUTINE) 

a) SOURCE 

b) REFLECTED 

c) DIFFRACTED FROM PLATE EDGE 

d) DIFFRACTED FROM FUSELAGE 
CHOPPED OFF EDGE 








PRINT, PLOT, AND/OR WRITE 
BINARY OUTPUT IN TERMS OF 
THETA AND PHI FIELD COMPONENTS 
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The subroutines for each of the scattered field components are all 
structured in the same basic way. First, the ray path is determined 
from the source to a particular scatterer and subsequently to the 
observation point using either the laws of reflection or diffraction. 
Each ray path, assuming one is possible, is then checked to see if it is 
shadowed by any structure along the complete ray path. If it is 
shadowed, the field is not computed and the code proceeds to the next 
scatterer or observation point. If the path is not interrupted, the 
scattered field is computed using the appropriate UTD solutions. The 
fields are then superimposed in the main program. This shadowing 
process is often speeded up by making various decisions based on bounds 
associated with the geometry of the structure. This type of knowledge 
is used wherever possible. 

The shadowing of rays is a very important part of the UTD 
scattering code. It is obvious that this approach leads to various 
discontinuities in the resulting pattern; however, the UTD diffraction 
coefficients are designed to smooth out the discontinuities in the field 
such that a continuous field is obtained. When a scattered field is not 
included in the result, the lack of its presence is apparent. This can 
be used to advantage in analyzing complicated problems. Obviously, in a 
complex problem not all the possible scattered fields can be included. 

In the UTD code the importance of the neglected terms are determined by 
the size of the so-called gliches or jumps in the pattern trace. 

If the gliches are small, no additional terms are needed for a good 
engineering solution. If the gliches are large, it may be necessary to 


10 



include more terms in the solution. In any case the user has a gauge 
with which he can examine the accuracy of the results and is not falsely 
led into believing a result is correct when in fact there could be an 
error associated with neglecting a higher order interaction term. 

The brief discussion of the operation of the scattering code given 
above should help the user get a feel for the overall code so he might 
better understand the code's capabilities and interpret its results. 

The code is designed, however, so that a general user can run the code 
without knowing all the details of its operation. Yet, he must become 
familiar with the input/output details which will be discussed in the 
next three sections. 

III. DEFINITION OF INPUT DATA 

The method used to input data into the computer code is presently 
based on a command word system. This is especially convenient when more 
than one problem is to be analyzed during a computer run. The code 
stores the previous input data such that one need only input that data 
which needs to be changed from the previous execution. Also, there is a 
default list of data so for any given problem the amount of data that 
needs to be input has been shortened. The organization of the input 
data is illustrated in Table II. 

In this system, all linear dimensions may be specified in either 
meters, inches, or feet; whereas, all angular dimensions are in degrees. 
All the dimensions are eventually referred to a fixed cartesian 


o 
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TABLE II 


BLOCK DIAGRAM OF THE INPUT DATA 
ORGANIZATION FOR THE COMPUTER CODE 
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coordinate system used as a common reference for the source and 
structure. The reference coordinate system is located at the origin of 
the composite ellipsoid. There is, however, a geometry definition 
coordinate system that may be defined using the "RT:" command. This 
command enables the user to rotate and translate the coordinates 
relative to the reference coordinate system and is used to input plate 
geometry in terms of the best coordinate system for this particular 
substructure. Once the "RT:" command is used, all the input following 
the command will be in that rotated and translated coordinate system 
until the “RT:" command is called again. The only exception to this is 
that the composite ellipsoid will always be in the reference coordinate 
system. See below for more details. There is also a separate 
coordinate system that can be used to define the pattern coordinates. 
This is discussed in more detail in Section III-C in terms of the "PD:" 
command . 

It is felt that the maximum usefulness of the computer code can be 
achieved using it on an interactive computer system. As a consequence, 
all input data are defined in free format such that the operator need 
only put commas or spaces between the various input variables. This 
allows the user on an interactive terminal to avoid the problems 
associated with typing in the field length associated with a fixed 
format. This method also is useful on batch processing computers. Note 
that all read statements are made on unit #5, i.e., READ (5,*), where 
the symbol refers to free format. Other machines, however, may have 
different symbols representing free format. 
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In all the following discussions associated with logical variables 
a "T" will imply true, and an “F" will imply false. The complete words 
true and false need not be input since most compilers just consider the 
first character in determining the state of the logical variable. 

The following list defines in detail each command word and the 
variables associated with them. Section VI will give specific examples 
using this input method. Note that the program halts execution by 
sensing the end-of-file mark associated with the input data stream. 
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COMMAND PART: 


A. Unit and Frequency Commands : 

Al. COMMAND UN: Set Linear Units Used for Input 

A2. COMMAND FQ: Frequency Input 

B. Fuselage Geometry Related Commands : 

Bl. COMMAND FG: Fuselage Geometry Input 

B2. COMMAND FB: Fuselage Blockage Modeled by Plates 

B3. COMMAND FC: Fuselage Chopped Off 

C. Source Geometry Related Commands : 

Cl. COMMAND SG: Source Geometry Input 

C2. COMMAND SP: Superposition Fields from Several Sources 

C3. COMMAND LS: Line Source Distribution Along Z-axis Used in 

Array Pattern 

D. Plate Geometry Related Commands : 

Dl. COMMAND PG: Plate Geometry Input 

D2. COMMAND PI: Initialize Number of Plates to be Retained 

E. Pattern Cut Related Commands : 

El. COMMAND PD: Conical Pattern Data Desired 

E2. COMMAND RT: Translate and/or Rotate Coordinates 

F. Specific Calculation Related Commands : 

FI. COMMAND TO: Test Data Generation Option 

F2. COMMAND RD: Reflection/Diffraction Included in Computation 

F3. COMMAND DD: Double Diffraction Included in Computation 

Reset Input Data to Default Case 


F4 


COMMAND RS 



G. Execute and Output Related Commands: 


Gl. 

COMMAND LP: 

Line Printer Listing of Results 

G2. 

COMMAND PP: 

Pen Plot of Results 

G3. 

COMMAND BO: 

Binary Outputs of E-THETA and E-PHI Pattern 
Results 

G4. 

COMMAND EX: 

Execute Program 
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. Unit and Frequency Commands: 


Al. COMMAND UN: 


READ: IUNIT 


This command enables the user to specify the units used for all 
following linear dimensions in the input data list. Note that this 
command should be defined before the other geometry -related commands. 


1. READ: I UN IT 


a) IUNIT: This is an integer variable that defines the units, 


1 -*• meters 


IUNIT = 


2 -*• feet 


3 > inches 


The default is 1. 
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A2. COMMAND FQ: 


© 

1 

/" 'READ: NFR'EQ.FREOI .DFREQ 1 

I 

© 


This command enables the user to specify the operating frequencies 
of the antennas. 

1. READ: NFREO, FREuI, DFREQ 

a) NFREQ: This is an integer variable that specifies the 

number of different frequencies. 


b) FREQI, DFREQ: They are real variables that specify the 

start and increment of the frequency loop, 
respectively, in Gigahertz. 

The default is N=l, FREQ I = . 2997925 , and DFREQ=. 2997925. 
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B. Fuselage Geometry Related Commands : 


Bl. COMMAND FG: 





/ READ: AX.BX.CX.DX / 


f 

/ READ: 

LSOUR / 

> 

f 

/ READ: (PVO(N), N=l,3) / 


I 



This command enables the user to model the fuselage by a composite 
ellipsoid. (See Section IV for more details in defining fuselage 
geometry.) 

1. READ: AX,BX,CX,DX 

a) AX,BX,CX,DX: These are real variables that specify the 

semi -mi nor axes of the composite ellipsoid 
used to model the fuselage as shown in Figure 
1. Note that the cross-section and profile 
of an ellipsoid is an ellipse. 
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2. READ: LSOUR 


a 


Figure 1 


) LSOUR: This is a logical variable defined by T or F. 

It is used to indicate if one wishes to 
define the location of the origin about which 
pattern is taken at the phase reference of 
antennas. 

The default is LSOUR = ‘FALSE* 



Definition of fuselage geometry. 
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3. READ: (PVO(N) , N-1,3) 


a) PVO(N): This is a real dimensioned array that defines the 

location of the origin about which the pattern is 
taken, i.e., PVO(N) = (x,y,z). Note PVO(N) is not 
used when LSOUR is .TRUE, but it must be input. 


B2. COMMAND FB: 



i 


/ READ: MPXFB / 


| DO 5201 MP=1 , MPXFB | 

i 

/ READ: MEXFB(MP) / 

1 

I DO 5201 ME=1 , MEFB | 

1 

/ READ: (PVFB(N,ME,MP) , N=l,3) / 

I 

5201 CONTINUE | 


I 



This command enables the user to model the fuselage blockage by 
plates. The blockage plates prevent the rays from passing through the 
fuselage and contributing to the field calculations. Note that these 
plates only shadow energy, i.e., they don't scatter energy. 
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1. READ: MPXFB 

MPXFB: This is an integer variable which defines the maximum 

number of plates to be used in modeling the fuselage 
blockage. MPXFB can not exceed 2! Normally, one uses 
two plates, i.e., one horizontal and the other 
vertical. Their contour is a best fit match to the top 
and side profiles of the fuselage. 

2. READ: MEXFB(MP) 

MEXFB(MP): This is a dimensioned integer variable which 

defines the maximum number of corners of each 
fuselage blockage plate. MEXFB(MP) can not exceed 
6 ! 

3. READ: (PVFB(N,ME,MP) , N=l,3) 

PVFR(N,ME,MP) : This is a triply dimensioned real 

variable. It is used to specify the 
location of the MEth corner of the MPth 
plate. It is input on a single line 
with the real numbers being the X,Y,Z 
coordinates of the corner which 
corresponds to N=l,2,3, respectively, in 
the array. For example, if the first 
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plate and 2nd comer is located at x=2, 
y=0, z=20, then it is represented by 


PVFB(1,2,1)=2. 

PVFB(2,2,1)=0. 

PVFB(3,2,1)=20. 

This data is input as: 2., 0., 20. 


B3. COMMAND FC: 


® 

1 

/ READ: LZC1.LZC2 / 

l 

/ READ: ZC1 ,ZC2 / 

i 

® 


This command enables the user to chop off the fuselage. This 
command is ueful in modeling the radome bulk/head portion of an aircraft 
fuselage. Using this command the fuselage ellipsoid is cut at right 
angles to the z-axis which forms an abrupt termination of the fuselage. 

1. READ: LZC1.LZC2 

LZC1/LZC2: These are logical variables defined by T or F. 

They are used to indicate if the fuselage will be 
chopped off in the ZC1/ZC2 location. 
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2. READ: ZC1, ZC2 


ZCI/ZC2: This is a real variable which defines positive/ 

negative Z location at which the fuselage is 
chopped, respectively. Note ZCI(ZC2) can be any 
number when LZC1(LZC2) is .FALSE. 
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This command enables the user to specify the location and type of 
antenna to be used. The geometry is illustrated in Figure 2. 

1. READ: PHS.ZS 

PHS,ZS: These are real variables used to specify the phi -angle 

(in degrees) and Z location of the antenna phase 
reference point. (Refer to Figure 2) Note: -90° < 

PHS < 90°. 

2. READ: MSX 

MSX: This is an integer variable which defines the maximum 

number of elemental radiators to be considered during 
execution of the program. Presently, 1 < MSX < 10. 

3. READ: RHOA(MS) , PHIA(MS) 

a) RHOA(MS): This is a dimensioned real variable which 

defines the distance that a single antenna 
element is positioned away from the antenna 
phase reference point. It is shown in Figure 3 
in terms of p/\. 

b) PHIA(MS) : This is a dimensioned real variable used to 

specify the angle (<t>/\ is in degrees) relative to 
the antenna coordinates shown in Figure 3. 
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ANTENNA PHASE 
REFERENCE POINT 






4. READ: SLOTAA(MS) , SLOTBA(MS), BETADA(MS) , SMONOA(MS) , 

JANTA(MS) 

a) SLOTAA(MS), SLOTBA(MS): These are real variables uesed to 

specify the narrow (parallel with E field) and 
broad (perpendicular to E field) dimensions of 
the slot in specified units. 

b) BETADA(MS): This is a real variable used to specify the 

angle (in degrees) of the slot relative to the 
fuselage axis. If BETADA=0, then it is an 
axial slot. If BETADA=90., then it is a 
ci rcuprfsrsnt i dl slot. 

c) SMONOA(MS): This is a real variable used to specify the 

length of the monopole in specified units. 

Note that SMONOA should not exceed a quarter 
wavelength. 

d) JANTA(MS): This is an integer variable used to specify 

the type of antenna considered in the 
computation: 

( 1 arbitrary oriented slot 
3 -► radial monopole. 
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5. READ: WMA(MS), WPA(MS) 


a) WMA(MS) , WPA(MS) : These are real variables used to specify 

the magnitude and phase (in degrees) of 
the excitation of the MSth antenna. If an 
array is used, then the excitation 
including the coupling effect on the 
radiators may be obtained using a 
thin-wire [4] as shown in the results 
section. 


C2. COMMAND SP: 



/ READ: LSUPER, MM, WP / 



This command enables the user to superimpose the fields calculated 
for several sources. But, one should note that this command can be used 
only when sources are operating at the same frequency. Note that the 
fields are superimposed with each execution of the program if LSUPER = 
TRUE. 

1. READ: LSUPER ,WM,WP 

a) LSUPER: This is a logical variable defined by T or F. It 

is to indicate if one wishes to superimpose fields 
or not. 
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b) WM,WP: These real variables are used to specify the 

magnitude and phase (in degrees) of the source 
relative to the first source in the superposition 
string. 

COMMAND LS: 



COMMAND LS (continued) 







COMMAND LS (continued) 




/ R£AD: SLENG, CTE 

W, NPOW, THSCAN / 


t 






/ READ: SLENG, N 

POW, RHOP, ZOP / 



This command enables the user to specify a line source distribution 
along the z-axis. It is used in an array pattern multiplication 
analysis. This command applies only when one has a uniform geometry 
along the axis of the fuselage. 

1. READ: LSTERM 

a) LSTERM: This is an integer variable that indicates the 

type of line source distribution treated. The 
current distribution and, therefore, the following 
inputs vary according to the following table. 
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LSTERM=0: I(z) - 


LSTERM=1 : 


LSTERM=2: 


LSTRM=3: 


LSTERM=4: 


READ: 

READ: 


NINP 

S |I N | e N «( z -(N-l/2)Az) 
N=1 

NINP.DELZ 

CURM.CPHAS 


I(z) 


r -otz -jB -ct 

= I e + re 


(2L-z) -jkz cos9 t 


READ : SLENG , DBATT ,GAMM .BETA , THSCAN 


2(z- 


LM] 


2 N 


I(z) = < ■ 

i 

READ: SLENG, NPOW, THSCAN 


-jkz cos0 


I (z) = 


^COS JTZ^ ^ + 


-jkz cos0 ( 


READ: SLENG, CTERM.NPOW.THSCAN 


I(z) = 


COS TTZ 

L 


-jk(p-Po) 


f 2 " 2 

where p = yp 0 + ( z_z o) 
READ: SLENG, NPOW.RHOP.ZOP 
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The input data is interpreted as follows: 


a) NINP: 


b) DELZ: 


c) CURM,CPHAS: 


d) SLENG: 


e) DBATT : 


f) GAMM.BETA: 


g) THSCAN: 


h) NPOW: 


This is an integer variable that defines the number 
of current samples. 

This is a real variable (Az) that defines the 
current sample spacing in wavelengths. 

These are real variables that define the magnitude 
( 1 1 IN) and phase (<|>N) of the current elements. 

This is a real variable (L) that defines the 
length of the linear array, in wavelengths. 

This is a real variable that defines the 
attenuation (in dB) along the total length (SLENG) 
of the array. 

Note that a is related to DBATT. 

These are real variables (r and s) that define the 
magnitude and phase (in degrees) of the reflection 
coefficient at the end of the traveling wave 
antenna (LSTERM=1). 

This is a real variable that defines the scan angle 
(in degrees) of the array. 

This is an integar variable (N) that defines the 
exponent in the previous equations. 
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i) CTERM: This is a real variable that defines the constant 

(C) in the previous equations. 

j) RH0P,Z0P: These are real variables that define the phase 

distribution across an aperture. Note that RHOP 
and ZOP are specified in wavelengths. In terms of 
the previous definition for the case (LSTERM=4) 
RH0P=p o and ZOP=Z 0 . 

D. Plate-Geometry Related Commands : 

01. COMMAND PG: 
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This command enables the user to define the geometry of the flat 
plate structures to be considered. The geometry is illustrated in 
Figure 4. It can be called repeatedly up to 25 times. If the higher 
order terms (R/R, R/D, D/R, D/R) are not defined, it can be extended up 
to 26. Note that each side of the chopped-off fuselage is automatically 
counted as three plates in the computer code when the "FC:" command is 
used. See Section IV for further details in defining the corner 
locations. 



Figure 4. Definition of flat plate geometry. 
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1. READ MCX(MPX) ,LATACH(MPX) 


a) MCX(MPX): This is a dimensioned integer variable. It 

is ued to define the number of corners (or 
edges) on the MPXth plate. Presently, 
1<MCX(MPX)<6 with 1<MPX<25. 

b) LATACH(MPX): This is a logical variable defined by T or F. 

It is used to indicate if the MPXth plate is 
attached to the fuselage (T) or not (F). If 
it is true, the first and last corners of 
attached plates should be specified on or 
near the fuselage. If they are not attached, 
the program will automatically attach the 
first and last corners. 


2. READ: (PVC(N,ME,MPX) ,N=1,3) 

As stated earlier, the locations of the corners of the flat plates 
are input in terms of the x, y, z coordinates in the specified cartesian 
coordinate system. 

a) PVC(N,ME,MPX) : This is a triply dimensioned real variable. 

It is used to specify the location of the 
MEth corner of the MPXth plate. It is input 
on a single line with the real numbers being 
the x, y, z coordinates of the corner which 
correspond to N=l,2,3, respectively, in the 
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array. For example, the array will contain 
the following for plate #1 and corner #2 
located at x=2, y=4 , z=6: 

PVC( 1,2,1) = 2 
PVC(2,2,1) = 4 
PVC(3,2,1) = 6 


This data is input as: 2., 4., 6. 

Considering the flat plate structure given in Figure 4, the input 
data is given by 


1 ., 1 ., 0 . 
- 1 ., 1 ., 0 . 
-l.,-l., 0. 

o. 


corner #1 
corner #2 
corner #3 
corner #4 


Presently, 1 _< MPX £ 25 
1 < ME <_ 6 
1 < N <3 

Note that the limits on the number of plates, and corners are only due 
to the size of the arrays. 



1 


/ Wm MPHOLD / 


D2. COMMAND PI: 



This command enables the user to specify those consecutive plates 
which will remain for the next calculation. Its useful when one 
simulates a complicated model by many plates and wants to know the 
effect of eliminating some plates from a configuration. The usage is 
illustrated in example 2. 

1. READ: MPHOLD 

MPHOLD: This is a real variable used to specify the number of 

plates to be retained for the next calculation. One 
should note that the first MPHOLD plates are retained. 
For example, if MPH0LD=2, then plates #1 and #2 remain 
in the input data list for the next computation. 

E. Pattern Cut Related Commands : 

El. COMMAND PD: 


© 

i 

/ READ: THC , PHC, THETA / 

V 


/ READ: IPS, IPF, IPP~7 


> 

f 

/ READ: i 

LFAR, R / 

1 

t 
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This command enables the user to define the pattern axis of 
rotation, the angular range, and the range from origin to receiver for 
the desired conical pattern. 

This set of data is associated with the conical pattern desired 
during execution of the program. The pattern axis is defined by the 
spherical angles (THC,PHC) as illustrated in Figure 5. These angles 
define a radial vector direction which points in the direction of the 
pattern axis of rotation. These angles actually set-up a new coordinate 
system in relation to the reference coordinates. The new cartesian 
coordinates defined by the subscript "p" are found by first rotating 
about the z-axis the angle PHC and, then, about the y n -axis the angle 

r 

THC. The pattern is, then, taken in the "p" coordinate system in terms 
of spherical angles. The theta angle of the pattern taken about the 
Zp-axis is defined by THETA. The phi angle is defined by the next read 
statement. In the present form the program will compute any conical 
pattern in which THETA is used as the conical pattern angle about the 
Zp-axis for the complete pattern calculation. 

As an aid in setting up the "p" coordinate system the following set 
of equations give the relationships between (THC, PHC) and the Xp, yp, 
and Zp-axes: 

x p = cos ( PHC ) s i n ( THC+90 0 ) x+s i n ( PHC ) s i n ( THC+90 0 )y+cos ( THC+90 0 ) z 

y p = cos ( PHC+90 0 ) x+s i n ( PHC+90® )y 
and 

Zp = cos(PHC)sin(THC)x+sin(PHC)sin(THC)y+cos(THC)z 
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where 0 £ THC <_ 180° and 0 _< PHC _< 360°. Note that the "p" axes are 
defined as radial vector directions in a spherical coordinate system. 

In its present form it should be noted that the user may not be able to 
define the Xp-axis at the starting location that he desired. In 
addition, the rotation of the pattern may be in the opposite sense using 
this approach. However, these problems can be easily overcome with 
properly written plot routines. 

1. READ: THC, PHC, THETA 

a) THC, PHC: These are real variables. They are input in 

degrees and define the axis of rotation about 
which a conical pattern will be computed (see 
Figure 5.). 

b) THETA: This is a real variable. It is input in degrees 

and ued to define the conical angle (9p) about 
the axis of rotation for the desired pattern. 

2. READ: IPS.IPF.IPD 

a) IPS,IPF,IPD: These are integer variables used to define 

angles in degrees. They are, respectively, the 
beginning, ending, and incremental values of the 
phi pattern angle (<|>p). 

As a result of the input given by the two previous read statements, 
the user has completely defined the desired conical pattern to be 
computed during execution of the program. 
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X 


Figure 5. Definition of pattern axis. 
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3. READ: LFAR.R 


a) LFAR: This is a logical variable defined by T or F. It is 

used to specify if the far field pattern is desired 
or not. 

b) R: This is a real variable which is used to define the 

range in linear units from the origin to the 
receiver. Note R can be any number when LFAR is 
.TRUE, in that it is not used in the calculation. 


E2. COMMAND RT: 





/'READ: (TR(N) > N=1 > 5T7 

yf 

/"READY 'THZR, PHZRT THXK~,~~RriXR / 



This command enables the user to translate and/or rotate the 
coordinate system used to define the input data in order to simplify the 
specification of the plate geometry. The geometry is illustrated in 
Figure 6. 
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1. READ: (TR(N) , N-1,3) 

a) TR(N): This is a dimensioned real variable. It is 

used to specify the origin of the new 
coordinate system to be used to input the data 
for the plate structures. It is input on a 
single line with the real numbers being the 
x,y,z coordinates of the new origin which 
corresponds to N=l,2,3, respectively. 

2. READ: THZR , PHZR, THXR, PHXR 

a) THZR, PHZR: These are real variables. They are input in 

degrees as spherical angles that define the 
ZR-axis of the new coordinate system as if it 
was a radial vector in the reference coordinate 
system (see Figure 6). 

b) THXR, PHXR: These are real variables. They are input in 

degrees as spherical angles that define the 
XR-axis of the new coordinate system as if it 
was a radial vector in the reference coordinate 
system (see Figure 6). 

The new x R -axis and z R -axis must be defined orthogonal to each 
other. If they are not, the program aborts with a warning. The new 
y R -axis is found from the cross product of the x R and z R axes. All the 
subsequent inputs will be made relative to this new coordinate system. 
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which is shown as (xr, yR, zr) unless command "RT:" is called again and 
redefined. It is always defined relative to the reference coordinate 
system. 

F. Specific Terms Related Commands : 

FI. COMMAND TO: 



/ READ: LDEBIlfi, LTEST, LOUT / 



True 


r READ: (LTRM(J), J=l,8) / 

I 

/ READ: MPI , MPF, MPS / 

1 

nxm~: — (MFKMp-JVWCMp-) 

/_ MES(MP), NP=MPI ,MPF ,MPS) 
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This command enables the user to obtain an extended output of 
various intermediate quantities in the computer code. This is useful in 


testing the program or in analyzing the contributions from various 
scattering mechanisms in terms of the total solution. 

1. READ: LDEBUG, LTEST, LOUT 

a) LDEBUG: This is a logical variable defined by T or F. It 

is ued to debug the program if errors are 
suspected within the program. If set true, the 
program prints out data on unit #6 associated with 
each of its internal operations. These data can, 
then, be compared with previous data which are 
known to be correct. It is, also, used to insure 
initial operation of the code. Only one pattern 
angle is considered, (normally set false) 

b) LTEST: This is a logical variable defined by T or F. It 

is used to test the input/output associated with 
each subroutine. The data written out on unit #6 
is associated with the window variables of the 
subroutine. They are written out each time the 
subroutine is called. It is, also, used to insure 
initial operation of the code. Only one pattern 
angle is considered, (normally set false) 
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c) LOUT: 


This is a logical variable defined by T or F. It 
is used to output data on unit #6 associated with 
the main program. It is also used to initially 
insure proper operation. It can be used to 
examine the various components of the pattern, 
(normally set false) 

2. READ: LTERM.LCORNR 

a) LTERM: This is a logical variable defined by T or F. It 

is used to tell the code whether or not individual 
terms are desired during the computation. 

(normally set false) 

b) LCORNR: This is a logical variable defined by T or F. It 

is used to tell the code whether or not corner 
diffraction is desired during the computation, 
(normally set true) 

3. READ: (LTRM(J), ,1-1,8) 

a) LTRM(J): These are logical variables defined by T or F to 
specify a set of individual scattering components 
that are to be included in the scattered field 
computation. The components are defined by the 
following number designations: 

J=l: source field 

J=2: single reflected field 
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J=3: single diffracted field 

J=4: diffracted field from chopped fuselage 

J=5: double reflected field 

J=6: reflected-diffracted field 

J=7: diffracted-reflected field 

J=8: double diffracted field 

The default values are LTERM(J) • 

(Note: To get the reflected-diffracted and/or double 

diffracted field one must accompany this command with COMMAND 
"RD:" and/or "DD:" , respectively.) 

4. READ: MPI,MPF,MPS 

a) MPI,MPF,MPS: These are integer variables to define the 

plates used in the computation, where 
MPI = initial plate 
MPF = final plate, and 

MPS = increment in plates going from initial 
to final plate, 

(Note: MP 1=1 , MPF=3, and MPS=2 imply plates 1 and 3 are 

included in the computation.) 

5. READ: (MEI(MP), MEF(MP), MES(MP), MP=MPI , MPF, MPS) 

a) MEI (MP) ,MEF(MP) ,MES(MP) : These are dimensioned integer 

variables to define the edges on the MPth 
plate used in the computation, where 
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MEI(MP) = initial edge on plate MP 
MEF(MP) = final edge on plate MP, and 
MES(MP) = increment in edges going from 
MEI(MP) to MEF(MP). 


F2. COMMAND RD; 



1 


/ READ: NRDX~7 

i 

/ DO 5622 NRD=1 ,NRDX / 

/ READ: MPIRD(NRD) , MP2RD(NRD) 7 

I 

/ 5622 CONTINUE"? 

I 



1. READ: NRDX 

NRDX: This is a real variable used to specify the number of 

reflection-diffraction terms desired. Presently, 

0 < NRDX < 40. 
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2. READ: MPIRD(NRD) , MP2RD(NRD) 

MPIRD(NRD): This is an integer dimensioned array used to 

specify the plate number from which the first reflection 
occurs. 

MP2RD(NRD): This is an integer dimensioned array used to 

specify the plate number from which the diffraction 
occurs. 

(Note: The usage of this command is illustrated in example 2.) 

F3. COMMAND DD: 



I 


/ READ: NDDX / 

I 

* /" DO 4022 NdP'=1,NDDx / 

i 

/ READ: MP1DD(NDD) > ME1DD(NDD).MP2DD(NDD) > ME2DD(NDD) / 

/ 4022 CONTINUE / 
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1. READ: NDDX 


a) NDDX: This is an integer variable that specifies the total 

number of double diffraction terms desired: 
Presently, 0 £ NDDX £ 10. 


2. READ: MPIDD(NDD) , MEIDD(NDD) , MP2DD(NDD) , ME2DD(NDD) 

a) MPIDD(NDD) , MEIDD(NDD) : These are integer dimensioned 

arrays used to specify the plate and edge number, 
respectively, from which the first diffraction 
occurs. 


u \ 
») 


unorvr» f \ ur onn/Mn\ . 


TU a *> M r\ n nn a m /ll mArtf 1 rtna/^ 

I I1C5C ai c i nuclei u i men j i wntu 


arrays used to specify the plate and edge number, 
respectively, from which a second diffraction 
occurs. 


(Note: The usage of this command is illustrated in example 2.) 


F4. COMMAND RS : 

This command enables the user to reset the input data to the 
default case. There is no input data associated with this command. 
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6. Execute and Output Related Commands : 
Gl. COMMAND LP: 




This command enables the user to obtain a line printer listing of 
the total fields (E 0pf E$p). 

1. READ: LWRITE 

LWRITE: This is a logical variable defined by T or F. It is 

used to indicate if a line printer output is desired 
or not. 


G2. COMMAND PP: 



/ READ: LPLOT 7 

/ READ : PLTNUM, RADIUS. TKT1 


This command enables the user to obtain a pen plot of the total 


fields (E 9p> E*p). 

1. READ: LPLOT 

LPLOT: This is a logical variable defined by T or F. It 

is used to indicate if pen plot is desired or 
not. 


2 . 


READ: PLTNUM, RADIUS, IPLT 

a) PLTNUM: This is a real variable used to indicate the type 

of polar plot desired, such that 


PLTNUM = 


l+E-THETA and E-PHI are plotted separately, 
j 2+E-THETA and E-PHI are plotted in the same 
pi ot. 

3+Both 1 and 2. 


b) RADIUS: This is a real variable that is used to specify 

the radius of the polar plot. 


c) IPLT: 


This is an integer variable that indicates the 
type of polar plot desired, such that 


IPLOT 


1 + field plot 
\ 2 power plot 

3 > dB plot 
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63. COMMAND BO: 



This command enables the user to obtain a binary output of the 
complex E-THETA and E-PHI patterns values. This output is useful if one 
wishes to input the aircraft code results into another program. This 
might be useful, for example, to study array patterns. One can run each 
array element individually, and then process the array pattern by 
appropriately adjusting the amplitude and phase of each element pattern. 
In this way numerous array patterns can be obtained without running the 
aircraft code excessively. 

1. READ: LBOUT 

a) LBOUT: This is a logical variable defined by T or F. It 

is used to indicate if the binary output is desired 
or not. The output format is specified within the 
source listing of the code. 


G4. COMMAND EX : 

This command is used to execute the code so that the total fields 
may be computed. After execution the code returns for another possible 
command word. 
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This concludes the definition of all the input parameters to the 
program. The program would, then, run the desired data and output the 
results on unit #6. However, as with any sophisticated program, the 
definition of the input data is not sufficient for one to fully 
understand the operation of the code. In order to overcome this 
difficulty the next section discusses how the input data is interpreted 
and used in the program. 


IV. INTERPRETATION OF INPUT DATA 


This computer code is written to require a minimum amount of user 


information such that the burden associated with a complex geometry w 


^11 
I \ \ 


be organized internal to the computer code. For example, the operator 
need not instruct the code that two plates are attached to form a convex 
or concave structure. The code flags this situation by recognizing that 
two plates have a common set of corners (i.e., a common edge). So if 
the operator wishes to attach two plates together he needs only define 
the two plates as though they were isolated. However, the two plates 
will have two identical corners. All the geometry information 
associated with plates having common edges is then generated by the 
code. The present code also will allow a plate to intersect another 
plate as shown in Figure 7(a). It is necessary that the corners 
defining the attachment be positioned a small amount through the plate 
surface to which it is being connected. Note that the edges of the two 
intersecting plates should be no closer than a quarter wavelength. 


59 



In defining the plate corners it is necessary to be aware of a 
subtlety associated with simulating convex or concave structures in 
which two or more plates are used in the computation. This problem 
results in that each plate has two sides. If the plates are used to 
simulate a closed or semi-closed structure, then possibly only one side 
of the plate will be illuminated by the antenna. Consequently, the 
operator must define the data in such a way that the code can infer 
which side of the plate is illuminated by the antenna. This is 
accomplished by defining the plate according to the right-hand rule. As 
one's fingers of the right hand follow the edges of the plate around in 
the order of their definition, his thumb should point toward the 
illuminated region above the plate. To illustrate this constraint 
associated with data format, let us consider the definition of a 
rectangular box. In this case, all the plates of the box must be 
specified such that they satisfy the right-hand rule with the thumb 
pointing outward as illustrated in Figure 7 (b). If this rule were not 
satisfied for a given plate, then the code would assume that the antenna 
is within the box as far as the scattering from that plate is 
concerned. 

In the "PG:" command, if LATACH(MPX)=T (i.e., the plate is attached 
to the fuselage), the program assumes that the first and last plate 
corners (PVC(N,1,MPX) and PVC(N,MCMX,MPX) ) are to be attached to the 
fuselage. The user must define the geometry accordingly. 

The plates can be attached to the ellipsoid as illustrated in 
Figure 7(c) and (d). However, when the plates are attached on the lower 
half of the ellipsoid, the y component of the first and last corners are 
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set equal to the y dimension of the ellipsoid center line as shown in 
Figure 7(e) and (f). It is important to note that the user need not 
exactly attach the first and last corners to the fuselage because the 
code will extend the edges and reset the first and final corner points 
on the fuselage as shown in Figure 7(g). 

In the "FG:" command, the composite ellipsoid is constructed from 
two ellipsoid sections positioned back to back and connected together 
such that its surface is continuous and smooth at the cross-section of 
the source location. The composite ellipsoid semi -major/mi nor axes are 
defined by AX, BX, CX and HX. The source location is defined by Z s and 
4> s . The case in which the source is positioned to the right of the 
coordinate system origin (Z s positive) as shown in Figure 8(a). It is 
assumed here that both the right and left ellipsoid coordinate systems 
are coincident. Then, the right side ellipsoid semi -major/mi nor axes 
and the source location are defined as 

(a F , b F , c F , V es , V rs ) = (AX, BX, CX, v es , V rs ) 

where 

AX sin (<f> s ) 

V es . • arcsln (Z 5 /CX) and V rs = Bx cos 
The parameters for the left side ellipsoid are given by the following: 

AX cos (VI 

i es 

a F = cos (v; s ) 
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1 BX cos ( v es) 


br- = 


cos (V es ) 


and 


C F - OX + Z $h 


where 


V „ = arcsin 
es 



CX cos v es +1 

1 

H 

1 


tan (V es ) "(DX + Z s ) 

M - 

1 


and 


'sh 


h - DX sin ( v es) 

I 

1 + sin (V es ) 


Note that Z s ^ is the distance between the right and left ellipsoid 
coordinate origins as shown in Figure 8. 

For the Case when the source is to the left of the origin (Z s 
negative as shown in Figure 8(b), the left ellipsoid semi -major/mi nor 
axes and source location are defined as 


(ap, bp, Cp, V e$ , V ps ) = (AX, BX, DX, V e$ , V ps ) 


where 

Z s AX sin (<|> s ) 

V es 58 arcsin M and v rs “ Wcos"T*V) 

The parameters for the right ellipsoid are given by 
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AX cos (V es ) 

I 

cos (V es ) 

BX cos (V es ) 

I 

cos (V es ) 


and 



where 


I 


V 


6 $ 



DX cos v es _1 

-1 


tan m + T7\ 

es ' - • “S' 

Mi +* 



and 


'sh 


Z s - CX sin (V es ) 
1 + sin ^ V es^ 
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Figure 7(a). Data format used to define a flat plate intersecting 
another flat plate. 



Figure 7(b). Data format used to define a box structure. 
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(e) 



Figure 7(c)-(f). Fuselage and wing geometries for aircraft model 

looking from the front. The antenna is assumed to be 
on the top portion of the models. 
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( A ) WHEN ZS > 0 > 


ft ft* 



co 


WHEN ZS < 0 , 



Figure 8. Composite ellipsoid geometry simulating the aircraft 
fuselage. 
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Using the "SG:" command, it is necessary that -90° £ PHS(MS) £ 

90°, In case the antenna is placed on the bottom part of the fuselage, 
the user must redefine the geometry such that <(> s PHS(MS) falls within 
the required angular range. This requires turning the aircraft 
upside-down. 

The code simulates fuselage blockage by using "FB:" COMMAND. If 
this command is activated, the code will determine if a ray strikes a 
fuselage blockage plate. If so, it will set that field component to 
zero. Thus, the shadowing effect of the fuselage can be simulated in 
this way. It is assumed that the higher-order diffraction and 
reflection fields from the fuselage are small in which case they are 
neglected. Thus, even though higher-order interactions between 
structures and the fuselage are not added in the computation, their 
absence will be apparent in the results. 

Finally, it must be kept in mind that the antenna should be kept at 
least a wavelength away from any diffracting edge. In fact all 
dimensions should be at least a wavelength. 

V. PROGRAM OUTPUT 

The basic output option from the computer code is a line printer 
listing of the results. If LWRITE=T in the input data list, the program 
will automatically generate a line printer output of the complex field 
values as shown in Figure 9. Recall that the results of the program are 
the E 0p and E^p radiation pattern values. In order to again describe 
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these pattern components, let us consider the various principal plane 
patterns treated in the previous section. The computer code allows for 
a rotation of coordinates such that one can take a pattern about the 
spherical angles (THC, PHC). The geometry that applies for each of the 
roll, elevation, and azimuth patterns used in the next section is 
illutrated in Figure in. Note that the 0p and <|>p angles are defined 
relative to the rotated pattern coordinates and that they change as THC 
and PHC are changed. Thus, E Q p is the theta component of the field 

4- A •> A 

(i.e., Egp=E • 9p) in the pattern coordinate system. Likewise, E^p=E • 4»p) • 
The total radiated electric field is denoted by £. 

In addition to the printed results, one has the option of obtaining 
a set of polar patterns. If LPLOT=T in the input data list, using the 
"PP:" command, the program will automatically plot the E Qp and polar 
patterns. These patterns are plotted such that the outer ring 
corresponds to the pattern maximum in each case. This polar plot 
routine was used to plot the data presented in the next section. 

One more output option is to get the binary output of the E^p and 
E^ patterns. If LBOUT =T in the input data list, using the "BO:" 
command, the program will automatically write the E 0 p and E^ results on 
unit number #11, i.e., WRITE (11). Note that this unit number increases 
one by one (i.e., #12, #13, ...) for each additional execution. This is 
a very useful output when one wishes to interface this program with 
another one. 
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Line printer output for the E 9p fields of Example 1. 
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Figure 9(b). Line printer output for the E^ p fields of Example 1. 
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(b) ELEVATION PLANE COORDINATES (THC*90 # , PHC« 90° 1 


Figure 10. Transformed coordinate systems for the conical pattern 
cuts. 


72 


A 2 p 



(c) AZIMUTH PLANE COORDINATES ( THC«90* PHC » 0° ) 


Figure 10. (Continued). 
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VI. APPLICATION OF CODE TO SEVERAL SIMPLE EXAMPLES 


The following two examples are used to illustrate some features and 
demonstrate the usage of the basic COMMANDS of the computer code. The 
effect of higher order terms in the solution is shown in example 2. 

Note that the patterns are plotted in decibels with each division being 
10 dB and that the labeling is not included. 

Example 1. Consider the radiation pattern of an antenna mounted on a 
composite ellipsoid for different pattern cuts. This 
example illustrates the usage of the COMMAND "FC:" and its 
effect on the pattern. The geometry is shown in Figure 11. 


x 




(b) TOP VIEW 

Figure 11. A monopole mounted on a composite ellipsoid. 
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The input data is given by: 


FG: 

5. . 6. .60. .20. 

F 

0. ,0. ,0* 

SG:MQNOFOLE 

25. . 3. 

1 

0 . , 0 . 

.4, .8,0., .25,3 

1.. 0. 

H): ROLL PLANE (FAR FIELD) 

0 . , 0 . , 80 . 

0,360,1 

T,1000. 

FP: 

T 

2 . 1 . . 3 

BC: 

PDrAZIMUm PLANE (FAR FIELD) 

80 . ,0 . , 80 . 

0,360,1 

T,1000. 

EX: 

PD: ELEVATION PLANE (FAR FIELD) 

80. ,80. ,80. 

0,360,1 

T,1000. 

EX: 

SG: MONOFQLE 

25., -10. 

1 

0 . , 0 . 

•4, .8,0., .25,3 

1.. 0, 

EX: 

FC: FUSELAGE CROPPED OFF 
F,T 

40., -14. 

ESC: 


The computed results are shown in Figure 12. 
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(c) 9 C =90°, <j> c =90° , 9 p =90° 



(d) 9 C =90°, 4> c =90° , 9 p =90° (e)9 c =90°, * c =90°, 9 p =90° 


Figure 12. Radiation pattern of monopole mounted on a composite 

ellipsoid at frequency .3 GHz. (a) (b) (c) source located 
at PHS=25°, ZS=3X (d) (e) source located at PHS=25°, 

ZS=10X and fuselage chopped off at ZC2=-14X for (e). 
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Example ?.: Consider the roll plane radiation pattern for a bent plate 


attached to a composite ellipsoid (5' x 6', 50' x 50'). The 
geometry is shown in Figure 13. The usage of "TO:" and 
"PI:" commands and their effect on the radiation pattern 
will be shown in this example. Various GTO terms involved 
in the computation are shown in Figure 14. 
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Figure 13. A bend plate attached to a composite ellipsoid. 


77 






The input data is given by: 


UN: FEET 
2 

FQ: 1GHz 
■ 1 , 1 ., 1 . 

F6: 

5.. 6..50..50. 

F 

OtfOafOt 

IB: 

2 

4 

4.5.0. .20. 

4.5.0. , -20. 

-4.5,0., -20. 

—4.5 ,0 . , 20. 

4 

0.,— 5.5,20. 

0.,5.5,20. 

0.,5.5,-20. 

0.,-5.5,-2Q. 

SG: MONORXE 
0 . , 0 . 

1 

0 . , 0 . 

.4, .0,0.,. 25,3 

1.. 0. 

PP: PEN PLOT 
T 

1,1.35,3 

PD: BOLL PLANE (NEAR FIELD) 
0.,0.,90. 

0,360,1 

F,1000. 

PG: 

4,T 

3. . 6. ,—20. 

3. . 9. ,—20. 

3. . 9 . .20. 

3. . 6 . . 20. 

PG: 

4,F 

3. . 9. ,-20. 

10.. 18. ,-20. 

10. ,18. ,20. 

3.,9.,20. 

BD: 

1 

1,2 


ED: 

4 

2.4. 2. 2 

1.4. 2. 2 

2.2. 1.4 

2 . 2 . 2. 4 

TO: TOTAL FIELD (INCLUDE DOUBLE TERMS) 

F,F,F 

T,T 

T,T,T,F,T,T,T,T 

1 , 2,1 

1,4,1 

1.4.1 
EX: 

TO: SOURCE FIELD ONLY 

F,F,F 

T,T 

T,F,F,F,F,F,F,F 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: REELECTED FIELD ONLY 

F,F,F 

T,T 

F,T,F,F,F,F,F,F 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: DIFFRACTED FIELD ONLY 

F,F,F 

T,T 

F,F,T,F,F,F,F,F 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: S+R 

F,F,F 

T,T 

T,T,F,F,F,F,F,F 

1 . 2.1 

1,4,1 

1,4,1 
EX: 
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TO: S+RtD (ONLY FIRST ORDER TERM INCLUDE) 

F,F,F 

T # T 

T,T,T,F,F,F,F,F 

1 , 2,1 

1 , 4,1 

1 . 4.1 
EX: 

TO: DOUBLE REFLECTION (R/R) 

F,F,F 

T,T 

F,F,F,F,T,F,F,F 

1 . 2.1 

1 , 4,1 

1 . 4.1 
EX: 

TO: S+RfR/R 

F,F,F 

T,T 

T,T,F,F,T,F,F,F 

1 . 2.1 

1 , 4,1 

1 . 4.1 
EX: 

TO: REFLEETION/DIFFRACTIGN TERM (R/D) 

F,F,F 

T,T 

F,F,F,F,F,T,F,F 

1 . 2.1 

1 , 4,1 

1 . 4.1 
EX: 

TO: R/RfR/D 

F,F,F 

T,T 

F,F,F,F,T,T,F,F 

1 . 2.1 

1 , 4,1 

1 . 4.1 
EX: 

TO: S+RHRR/R+R/D 

F,F,F 

T,T 

T,T,T,F,T,T,F,F 

1 . 2.1 

1 , 4,1 

1 , 4,1 
EX: 
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TO: DIFFRACTICN/REILECTIGN TEPH (D/P) 

F,F,F 

T,T 

F,F,F,F,F,F,T,F 

1 , 2,1 

1,4,1 

1.4.1 
EX: 

TO: S+IHCH-R/ Rf R/EH-D/R 

F,F f F 

T,T 

T,T,T,F,T,T,T,F 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: DOUBLE DIFFRACTION TERM (D/D) 

F,F,F 

T,T 

F,F,F,F,F,F,F,T 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: D/R+D/D 

F,F,F 

T,T 

F,F,F,F,F,F,T,T 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: ALL DOUBLE TERM 

F,F,F 

T,T 

F,F,F,F,T f T,T,T 

1 . 2.1 

1,4,1 

1.4.1 
EX: 

TO: 

F,F,F 

T,T 

T,T,T,F,F,F,F,F 

1 . 2.1 

1,4,1 

1,4,1 


FI: TAKE OFF SEOCND ELATE 
1 

FG: ADD ONE ELATE 
4,T 

20i 

3., -9. ,20. 

3., -9. f -20. 

3., -6.,-20. 

PP: 

T 

1.. 2..3 
EX: 
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The computed results are shown in Figures 15 and 16. Note that 
each pattern in Figure 15 is normalized to the same level so that one 
can see the relative significance of each term. An interesting result 
is shown in Figure 15(e) where the source and the reflected field are 
superimposed. These two terms form the classical "Geometrical Optics" 
(GO) solution. However, one should note that the GO solution is far 
from being complete as one can observe from the discontinuities in the 
pattern. Even when the first order terms of the GTD solution are 
superimposed, as shown in Figure 15(f), the pattern is still rough. 
Therefore, higher order terms which are in Figures 15(g) to 15(p) are 
added to eliminate the discontinuities. The final result is shown in 
Figure 15(a). It is clear that these higher terms can be significant in 
certain regions of the pattern as shown in Figure 15(p). In the last 
execution the geometry is modified by removing one plate and adding 
another as shown in Figure 16. 
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(a) total solution 

( S+R+D+R/R+R/D+D/R+D/D ) 



(c) reflected field (R) 


Figure 15. Roll plane radiation 


(b) source field (S) 



(d) diffracted field (D) 


pattern. 
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(e) S + R 


(f) S + R + D 



(g) reflected/reflected (h) S + R + R/R 

field (R/R) 

Figure 15. (Continued). 
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c-B- 


(f) reflected/diffracted 
field (R/D) 



(k) S + R + D + R/R + R/D 


(j) R/R + R/D 



(1 ) di f f racted/ref 1 ected 
field (D/R) 


Figure 15. (Continued) 



(m) S + R + R/R + R/D + D/R 



(o) D/R + D/D 



(n) diffracted/diffracted 
field (D/D) 



(p) second order interaction 
GTD terms 

(R/R + R/D + D/R + D/D) 


Figure 15. (Continued). 
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0 >) 

Figure 16. Total solution (S+R+D) after using "PI:" and "PG:“ 
commands. 
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VII. APPLICATION OF CODE TO AIRCRAFT SIMULATIONS 


To begin any simulation of an aircraft, one needs to start with a 
set of scale model drawings. A typical aircraft model consists of a 
composite ellipsoid fuselage plus flat plates simulating the other 
structures such as wings, stabilizers, etc. One can also use the 
"COMMAND FC:" to model the radome part of the aircraft. The radome is 
constructed of low dielectric constant material such that it is assumed 
to be totally transparent in these calculations. 

A wide variety of aircraft computer models such as commercial, 
private, military aircrafts and the space shuttle are given in the 


roi lowing examples. They serve to illustrate 




Capability, and 


validity of this general numerical solution. 

To begin the simulation procedure, one first finds the composite 
ellipsoid parameters for the aircraft fuselage. The ellipsoid surface 
should simulate the fuselage surface as accurately as possible near the 
antenna location. Once the composite ellipsoid dimensions are specified 
the plates are added to the model. 

This code allows for two different methods for defining one plate 
to be attached to another: 1) edge to edge attachment and 2) edge to 

surface attachment. Edge to edge attachment, as illustrated in Figure 
17 often requires that a plate edge be defined as two or three col inear 
edges as the program identifies this mode of attachment only by finding 
two identical pairs of corners. Note that the corners must be 
consecutive on both plates which means there actually exists an edge 
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Figure 17. Edge to edge plate attachment. 

between them. In the case to surface attachment, one plate is defined 
as penetrating a short distance through the surface of the second plate 
as illustrated in Figure 18. The program then defines the new junction 
edge and eliminates the smaller portion of plate #1 behind plate #2. 
Here care must be taken to assure that the new junction edge is 
completely contained within the bounds of plate #2, and no where nearer 
than a quarter wavelength or so to an edge of plate #2. 



Figure 18. Edge to surface plate attachment. 
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One thing which should be noted is that the attaching of the wings, 
stabilizers, and plates to the fuselage is automatically done by the 
computer as illustrated in Figure 7(g), Thus, the user need not worry 
about generating the correct input data to perform this task. 
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Example 3: Simulation of Boeing 737 


In this example, monopole and slot antennas are mounted on the 
fuselage of a Boeing 737 aircraft at various stations as shown in Figure 
19. 


A A/4 monopole mounted at station 220 just above the cockpit as 
shown in Figure 19. The line drawing of the 737 is shown in Figure 20 
and the computer model based on the input data is shown in Figure 21. 
The input data is as follows: 


UN: INCHES 
3 

FQ: 3.18 GHZ 
1,3.18,1. 

FG: BOEING 737 (STATION 220) 

77. . 74. .830. .308 .56 
F 

o.,o.,o. 

SG: NONOPOLE 
0.,-278. 

1 

0 . ,0 . 

1.537. 3 .074.0. .. 928525. 3 

1 .. 0 . 

PG: RIGHT WING 

4,T 

1.. 75.. 67 .952 

1.. 536.93.316. 14 

1.. 536. 93 .379.86 

1.. 75.. 240 .26 
PG: LEFT WING 
4,T 

1., -75., 240 .26 

1., -536.93,379.86 

1., -536.93,316. 14 

1., -75., 67. 952 

PG: VERTICAL STABILIZER 
4,T 

77. . 8.25. 618.55 
284.147,8.25,819.056 

284. 147.0. . 683 .696 

77. . 0.. 483. 19 


PG: VERTICAL STABILIZER 
4,T 

77. . 0.. 483. 19 

284. 147.0. . 683 .696 
284.147,-8.25,819.056 

77. , -8.25,618. 55 
PG: NOSE 

4.T 

5. 5, -10., -308. 56 
-.1,-10. ,-321 .6 
-.1,0. ,-321 .6, 

5.5.0. ,-308.56 
PG: NOSE 

4,T 

5. 5.0. , -308. 56 
-.1,0. ,-321 .6 
-.1,10. ,-321 .6 

5.5.10. , -308 .56 

PP: POLAR PLOT IN DB 
T 

U. 5,3 

PD: ROLL PLANE 
0.,0.,89.8 
0,360,1 
F,6000. 

EX: 

PD: ELEVATION PLANE 

90. . 90. .89.8 
0,360,1 
F,6000. 

EX: 
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PD: AZIMUTH PLANE 

PD: AZIMUTH PLANE 

93., 0., 89. 6 

93 . ,0 . ,80 . 

0,360,1 

0,360,1 

F,6000. 

F,6000. 

EX: 

EX: 

PD: AZIMUTH PLANE 

PD: AZIMUTH PLANE 

93.,0.,30. 

93. ,0., 100. 

0,360,1 

0,360,1 

F,6000. 

F,6000. 

EX: 

EX: 

PD: AZIMUTH PLANE 

PD: AZIMUTH PLANE 

93 . ,0 . ,60 . 

93. ,0., 110. 

0,360,1 

0,360,1 

F,6000. 

F,6000. 

EX: 

EX: 

PD: AZIMUTH PLANE 

PD: AZIMUTH PLANE 

93.,0.,70. 

93. ,0., 120. 

0,360,1 

0,360,1 

F,6000. 

F,6000. 

EX: 

EX: 


The three principle plane results are shown in Figures 22 to 24 and 
found to be in very good agreement with measurements. The experimental 
work was performed by the technical staff at NASA (Hampton, Virginia) 
using a 1/ 1 1th scale model of a Boeing 737 aircraft. It is noted that 
the measured results have some asymmetry in the patterns. This could be 
attributed to misalignment of the monopole with respect to the surface 
normal or the movement of the model due to shifting weight during the 
measurement. This misalignment of the monopole (approximately 3° tilted 
to the nose section from the actual surface normal) was detected by the 
calculated evaluation plane pattern and various azimuth plane patterns. 
To compensate this misalignment in the comparison with measured 
patterns, the conical pattern axis rotated 3° to the nose section was 
used in the following pattern calculations for the Station 220 case. 

The various azimuth plane patterns for this antenna location are 
computed and shown in Figures 25(a) to (g). In each case, the 
calculated results compare very favorably with the measurements. 
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ANTENNA LOCATION 

A - STATION 220 
B- STATION 250 
C - STATION 305 



Figure 19. Test locations for the antenna installation on the Boeing 
737 aircraft. 
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(c) TOP VIEW 


Figure 20. Boeing 737 aircraft. 
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(a) SIDE VIEW 



Figure 21. Computer simulated model of a Boeing 737 aircraft. The 
antenna is located at Station 220. 
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CALCULATED 

MEASURED 


TOP 



(SCALE: EACH 0 I V I S I ON- 1 OOB) 


Figure 22. Roll plane pattern of a X/4 monopole mounted at Station 220 
on top of a Boeing 737 aircraft. 
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TOP 



BOTTOM 

ISCRLEs EACH 0 I V I S I ON- 1 OOB) 


Figure 23. Elevation plane pattern of a A/4 monopole mounted at Station 
220 on top of a Boeing 737 aircraft. 
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(SCALE < EACH DIVISION- 10DB) 


Figure 24. Azimuth plane pattern of a X/4 monopole mounted at Station 
220 on top of a Boeing 737 aircraft. 
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NOSE 


LEFT 

WING 



- CALCULATED 

— MEASURED 


RIGHT 

MING 


TAIL 

(a) 0 p *5O° 


ISCALE, EACH OIVISION-IODB) 
NOSE 


LEFT 

HING 



RIGHT 

HING 


tail 

(b) 0p«6O° 


Figure 25. Azimuthal conical patterns of a X/4 monopole mounted at 
Station 220 on top of a Boeing 737 aircraft. 
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NOSE 



TAIL 

(d) 0 p «80° 


Figure 25. (Continued). 
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NOSE 



TAIL 


(f) 0 p «IIO* 

Figure 25. (Continued). 
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CALCULATED 

MEASURED 


NOSE 



TAIL 

(SCALE* EACH DIVISION- 100B) 

(g) 0p ■ 120° 


Figure 25. (Continued). 


103 


The next two cases, a X/4 monopole is mounted at stations 250 and 
305 on top of the fuselage of the Boeing 737 aircraft. The input data 
is as follows: 


UN: INCHES 
3 

FQ: 3.18 GHZ 
1.3.18,1. 

FG: BOEING 737 

77. . 74. .830. .308 .56 
F 

0 . , 0 . , 0 • 

SG: MONOPOLE (STATION 250) 
2.9,-248. 

1 

0 . ,0 . 

1 .537.3 .074.0. . .928525.3 

1 .. 0 . 

PG: VERTICAL STABILIZER 
4,T 

77. . 8.25.618.55 
284.147,8.25,819.056 

284. 147.0. .683 .696 

77. . 0.. 483. 19 

PG: VERTICAL STABILIZER 
4»T 

77. . 0.. 483. 19 

284. 147.0. . 683 .696 
284.147,-8.25,819.056 

77. , -8.25,61 8.55 


PG: NOSE 
4,T 

5 .5,-10. ,-308. 56 
-.1,-10. ,-321 .6 
-.1,0. ,-321 .6, 

5.5,0. ,-308 .56 
PG: NOSE 
4,T 

5. 5.0. ,-308. 56 
-.1,0. ,-321 .6 
-.1,10. ,-321 .6 

5.5. 10. , -308 .56 

PP: POLAR PLOT IN OB 
T 

1 ,1 .5,3 

PD: ELEVATION PLANE 

90. . 90. .90, 

0,360,1 

F,6000. 

EX: 


The only difference in the input data for the case when the 
monopole is at station 305 is in the specification of the source 
geometry. The source geometry for station 305 is given by the 
following: 


SG: MONOPOLE 737 (STATION 305) 

0. ,-193. 

1 

0 . , 0 . 

1.537.3.074.0. ..928525.3 

1.. 0. 
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It is noted that the antenna at Station 250 is mounted 4" off the 
fuselage centerline. Both calculated and experimental results for 
Stations 250 and 305 are presented in Figures 26 and 27, respectively. 
The results reveal good agreement between the theoretical predictions 
and scale model measurements. 
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BOTTOM 

(SCRLEs EACH OIVISION-IODB) 


Figure 26. Elevation plane pattern of a X/4 monopole mounted at Station 
250 (off center) on top of a Boeing 737 aircraft. 
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NOSE 


BOTTOM 

(SCALE: EACH OIVISION-IOOB) 


Figure 27. Elevation plane pattern of a X/4 monopole mounted at Station 
305 on top of a Boeing 737 aircraft. 
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The next antenna location considered on the Boeing 737 aircraft, is 
station 222 on the bottom of the fuselage just behind the nose. The 
computer simulated model for this case is shown in Figure 28. The input 


data is as follows: 


UN: INCHES 
3 

FQ: 3.18 GHZ 
1 ,3.18,1 . 

FG: BOEING 737 (STATION 222) 

66. . 55. .785. .232. 52 
F 

o.,o.,o. 

SG: MONOPOLE 
0., -144.6 
1 

0 . , 0 . 

1.537 .3 .074.0. ..928525.3 

1 .. 0 . 


FC: FUSELAGE CHOPPED-OFF 
F,T 

0., -186.97 

PP; POLAR PLOT IN DB 
T 

1,2.46,3 

PD: ELEVATION PLANE 

90. ,90. ,90. 

0,360,1 

F,6000. 

EX: 


Axial and circumferential slot antennas are also analyzed for 
station 222. The source geometry input data for these two cases are 
given by the following: 


SG: AXIAL SLOT 
0.,— 144.6 
1 

0 . , 0 . 

1.537.3. 074.0. ..928525.1 

1 .. 0 . 


The calculated elevation plane patterns for the monopole, axial 
slot and circumferential slot at station 222 are compared with measured 
patterns in Figures 29(a) through (c), and they are found to be in very 
good agreement. 


SG: CIRCUMFERENTIAL SLOT 
0., -144.6 
1 

0 . , 0 . 

1 .537 ,3 .074,90. ,.928525,1 

1 ., 0 . 
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RADOME 



Figure 28. Computer simulated model for a X/4 monopole mounted at 

Station 222 on the bottom of the fuselage of a Boeing 737 
ai rcraft. 
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CALCULATED 



BOTTOM 

(SCALE i EACH OIVISION-lOOBl 


(a) a x/4 monopie case 


Figure 29. Elevation plane patterns of an antenna mounted at Station 
222 on the bottom of a Boeing 737 aircraft. 
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CALCULATED 
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(SCALE* EACH DIVISION-IODB) 


(b) an axial slot case 


Figure 29. (Continued). 
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calculated 


TOP — --MEASURED 



BOTTOM 

ISCRLEi EACH OIVISION-IODB) 


(c) a circumferential slot case 


Figure 29. (Continued). 
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Finally, a x/4 monopole is located at station 950 on the bottom of 
the fuselage at the rear of the Boeing 737 aircraft. The computer 
simulated model for this case is shown in Figure 30. The input data is 
as follows: 


UN: INCHES 
3 

FQ: 3.18 GHZ 

1.3.18.1 . 

FG: BOEING 737 (STATION 950) 

77.. 74.. 580.. 200. 

F 

o.,o«,o. 

SG: MONOPOLE 
0.,300. 

1 

o.,o. 

1 .537 ,3. 074,0.,. 928525,3 

1 .. 0 . 

PG: RIGHT WING 
4,T 

8. 1.7.5. , -118. 
33.9,536.93,140. 

40.. 536. 93. 201. 

30.6. 75. . 107. 

PG: LEFT WING 
4,T 

30.6,-75. ,107. 

40., -536. 93, 201. 
33.9,-536.93,140. 

8. 1, -75. ,-118. 


PG: HORIZONTAL STABILIZER 

4*T 

-18.1 ,66., 443. 

-6 .3 ,207., 561. 

-0.4, 207., 620. 

—5. ,66. ,574. 

PG: HORIZONTAL STABILIZER 

4,T 

-5., -66., 574. 

-0.4,-207. ,620. 

-6 .3, -207., 561. 

-18.1 ,-66., 443. 

PP: POLAR PLOT IN DB 
T 

1,2.38,3 

PD: ELEVATION PLANE 

90. ,90. ,90. 

0,360,1 

F,6000. 

EX: 


The elevation plane radiation pattern is computed and found to be 
in very good Agreement with the measured pattern as shown in Figure 31. 
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ANTENNA LOCATION 


Figure 30. Computer simulated model for a X/4 monopole mounted at 

Station 950 on the bottom of the fuselage of a Boeing 737 
ai rcraft. 




BOTTOM 

(SCALE t EACH DIV ISION-IOOB) 


Figure 31. Elevation plane pattern of a x/4 monopole mounted at Station 
950 on the bottom of a Boeing 737 aircraft. 
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Example 4: Simulation of the KC-135 Aircraft 

In this example a monopole, axial slot and circumferential slot 
antennas mounted both over and forward of the wings are studied on a 
KC-135 aircraft. The line drawings of the KC-135 are shown in Figure 
32. The computer simulated models based on the input data are shown in 
Figures 33 and 34. 

The input data for the monopole mounted over the wings is as 


follows: 

UN: INCHES 

PG: VERTICAL STABILIZER 

3 

4,T 

FQ: 34.92 GHZ 

2. 946,0., 49. 492 

1 ,34.92,1 . 

14.076,0. ,58.025 

FG: KC-135 FUSELAGE GEOMETRY 

14. 076,-. 5, 64 .205 

3.3 ,3. ,72. ,8. 

2. 946 ,-.5, 55. 672 

F 

PG: NOSE 

0.,0.,18.81 

4,T 

SG: MONOPOLE OVER WING 

1. 39,-1 .36,-7.35 

0., 18.81 

1 .275,-1 .36,-9. 

1 

1.37,0., -9. 

o.,o. 

1.485,0., -7 .35 

.140,.2C0,90.,.0845,3 

PG: NOSE 

1.,0. 

4>T 

PG: RIGHT WING 

1.485,0., -7 .35 

4.T 

1.37,0., -9. 

-.5,3. ,12.31 

1.275,1.36,-9. 

-.5,28.5,36.41 

1.39,1.36,-7.35 

-.5,28.5,40.41 

PO: ELEVATION PLANE 

-.5,3. ,24.61 

90. ,90. ,90. 

PG: LEFT WING 

0,360,1 

4,T 

T,1000. 

-.5,-3. ,24. 61 

PP: PEN PLOT 

-.5,-28.5,40.41 

T 

-.5,-28.5,36.41 

1,1.71,3 

-.5,-3. ,12.31 

EX: 

PG: VERTICAL STABILIZER 

PD: ROLL PLANE 

4,T 

0 .,0. ,90 . 

2. 946,. 5, 55. 672 

0,360,1 

14.076,. 5, 64.205 

T,1000. 

14.076,0. ,58. 025 

EX: 

2.946,0. ,49.492 

PO: AZIMUTH PLANE 

90. ,0., 90 .2 

0,360,1 

T,1000. 

EX: 
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(c) TOP VIEW 


)del of a KC-135 aircraft. The antenna 
' the wings. 


9 


The radiation patterns of different antenna configurations on the 
KC-135 may be obtained simply by changing the source geometry command. 
The other source locations are specified as follows: 


SGt MONOPOLE FORWARD OF WINGS 
0 ., 8.34 
1 

o.,o. 

. 140 .. 280 . 90.,. 0845 ,3 

1 ., 0 . 


SGt AXIAL SLOT FORWARD OF WINGS 
0 ., 8.34 
1 

0 ., 0 . 

. 140 . . 280 . 0 . ,. 0843,1 

1 ., 0 . 


SGt CIRCUMFERENTIAL SLOT FORWARD 
0 ., 8.34 
1 

0 . ,0 . 

. 140 , . 280 . 90 . ,. 0843,1 

1 ., 0 . 


SGt AXIAL SLOT OVER WINGS 
0 ., 18.81 
1 

o.,o. 

• 140 , . 280 , 0 . ,. 0843,1 

1 ., 0 . 


SGt CIRCUMFERENTIAL SLOT OVER WINGS 
0 ., 18.81 
1 

0 . , 0 . 

. 140 ,. 280 , 90 . ,. 0843,1 

1 ., 0 . 


WINGS 


The elevation, roll and azimuth plane patterns for a short 
monopole, a circumferential KA-band waveguide and an axial KA-band 
waveguide mounted forward and over the wings are shown in Figures 35 to 
43. 

The computed results are found to be in very good agreement with 
the measurements in the elevation and roll planes. The precision 
pattern measurements (elevation and roll plane patterns) using the 1/25 
scale model were taken at NASA (Hampton, Virginia). Measured data was 
not available for the azimuth plane. 
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(a) Antenna mounted forward of wings 


TOP 



BOTTOM 

ISCRLEt EACH DIVISI0N-100B) 


(b) Antenna mounted over wings 


Figure 35: Elevation plane pattern for a X/4 monopole mounted on a 

KC-135 aircraft. 
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Figure 36. Elevation plane pattern for a circumferential KA-band 
waveguide mounted on a KC-135 aircraft. 
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(b) Antenna mounted over wings 


Figure 37. Elevation plane pattern for an axial KA-band waveguide 
mounted on a KC-135 aircraft. 
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Figure 38. Roll plane pattern for a X/4 monopole mounted on a KC-135 
ai rcraft. 
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Figure 39. Roll plane pattern for a KA-band circumferential waveguide 
mounted on a KC-135 aircraft. 
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Figure 40 
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. Roll plane pattern for a KA-band axial waveguide mounted on 
a KC -135 aircraft. 
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Figure 41. Azimuth plane pattern for a A/4 monopole mounted on a KC-135 
ai rcraft. 
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(a) Antanna aountad forward of wings 
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(b) Antanna aountad ovar wings 


. Azimuth plane pattern for a KA-band circumferential 
waveguide mounted on a KC-135 aircraft. 
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Figure 43. Azimuth plane pattern for a KA-band axial waveguide mounted 
on a KC-135 ai rcraft. 
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Example 5: Simulation of a Lindberg Antenna Mounted on a KC-135 

A Lindberg crossed slot antenna mounted on the fuselage of a KC-135 
is studied in this example. The line drawing of the KC-135 aircraft was 
shown in Figure 32. The computer simulated model based on the input 
data is shown in Figure 44. The input data is as follows: 


UN; INCHES 
3 

FQ: 6.25 GHZ 
1.6.25,1. 

F6: KC-135 FUSELAGE GEOMETRY 

3.3 .3. . 72. .8. 

F 

0.,0., 18.81 

SG: LINDBERG CROSSED-SLOT 

0. .2.25 
2 

Q.,0. 

.07375, 1.475.0. ,.0845,1 

1 .. 0 . 

0 . ,0 . 

.07375,1 .475, 90.,. 0845,1 

1 . 90; 

PG: RIGHT MING 
4,T 

-.5 ,3. ,12.31 
-.5,28.5,36.41 
-.5,28.5,40.41 
-.5 ,3. ,24.61 
PG; LEFT MING 
4.T 

-.5,-3. ,24.61 

-.5.-28.5,40.41 

-.5,-28.5,36.41 

-.5 ,-3., 12.31 

PG; VERTICAL STABILIZER 

4,T 

2. 946.. 5. 55. 672 

14.076. . 5.64 .205 

14.076.0. .58.025 

2.946.0. .49.492 


PG; VERTICAL STABILIZER 

4,T 

2.946.0. .49.492 

14.076.0. .58.025 
14. 076,-. 5, 64 .205 
2.946.-. 5, 55 .672 
PG: NOSE 

4,T 

1. 39,-1 .36,-7.35 
1.275.-1 .36,-9. 

1.37.0. , -9. 

1.485.0. ,-7.35 
PG: NOSE 

4.T 

1.485.0. . -7 .35 

1.37.0. . -9. 

1.275.1.36, -9. 

1.39.1.36, -7.35 
PC: ELEVATION PLANE 

90. . 90. .90. 

0,360,1 

T,1000. 

PP: PEN PLOT 
T 

1 .1 .625,3 
EX: 

PD: ROLL PLANE 
0.,0.,90.2 
0,360,1 
T,1000. 

EX: 

PD; AZIMUTH PLANE 

90.. 0..90. 

0,360,1 

T.1000. 

EX; 

PD; 45°C0NICAL CUT 

90.. 0..45. 

0.360.1 

T,1000. 

EXt 
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Various calculated patterns along with the measured results taken 
from reference [17J are presented in Figures 45 to 48. Again, good 
agreement is obtained. The gain level in each case is adjusted to 
compare with measurements. The E 9 pattern corresponds to the vertical 
component, to the horizontal component and E C p to the circularly 
polarized field. 
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Figure 45. Elevation plane pattern for Lindberg antenna mounted on a 
KC-135 aircraft. 
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Figure 45. (Continued). 
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Figure 46. Roll plane pattern for Lindberg antenna mounted on a 
KC-135 aircraft. 
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Figure 46. (Continued). 
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Figure 47. Azimuth plane pattern for Lindberg antenna mounted on a 
KC-135 ai rcraft. 


137 


NOSE 



TAIL 

(SCALE: EACH DIVISION-IOOB) 


M hp 


Figure 47. (Continued). 
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Figure 48: Azimuth conical pattern 0 p =45°) for Lindbert antenna 

mounted on a KC-135 aircraft. 
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Figure 48. (Continued). 
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Example 6: Simulation of Cessna 402B 


In this example a monopole antenna is mounted on the fuselage of a 
Cessna 402B aircraft just above the cockpit. The line drawing of the 
Cessna 402R is shown in Figure 49, and the computer model based on the 
input data is shown in Figure 50. The input data is as follows: 


UN: INCHES 
3 

FQ: FIVE GIGAHERTZ 

1.5.. 1. 

FG: CESSNA 402B 
8. 2 ,26. ,285. ,152. 

F 

0 . ,0 . ,0 . 

SG: MjNOEQLE 

0.,-10. 

1 

0 . , 0 . 

.414, .828,0. ,.25,3 

1.. 0, 

FG: LEFT WING (INNER EftRT) 
4,T 

-41.8,26.0,0. 

— 36 . 8 , 68 . , 0 . 

—36 .8 ,68 . ,65 . 

-41 .8 ,26., 65. 

FG: RIGHT WING (INNER EftRT) 
4,T 

-41 .8 ,-26., 65. 

—36 .8 ,-68. ,65. 

—36 .8 ,—68 . ,0 . 

-41.8,-26.0,0. 

FG: LEFT ENGINE (RIGHT SHE) 
6,F 

— 36 . 8 , 68 . , 0 . 

-36.8 ,68., -40. 

-21.8 ,68., -40. 

—21 .8 ,6 8 . ,85 « 

-36 .8 ,68., 85. 

-36 .8 ,68., 65. 


FG OUGHT ENGINE (USFT SHE) 

«,F 

—36 .8 ,-68. ,0 . 

-36 .8 ,-68., 65. 

-36 .8,-68. ,85. 

-21 .8 ,-68. ,85. 

-21 .8 ,-68. ,-40 • 

-36 .8 ,-68. ,-40 . 

FGi LEFT ENGINE (TOP) 

4,F 

-21.8 ,68., -40. 

-21 .8 ,106., -40. 

-21 .8 ,106 . ,85 . 

-21 .8 ,68., 85. 

FG OUGHT ENGINE (TOP) 

4,F 

-21 .8 ,-68., -40. 

—21 .8 ,—68 . ,85 . 

-21 .8 ,-106., 85. 

-21 .8 ,-106., -40. 

FG: LEFT ENGINE (LEFT SHE) 
6 F 

-36 .8 ,106., 65. 

-36 .8 ,106., 85. 

-21 .8 ,106 . ,85 . 

-21 .8 ,106 . ,-40 . 

-36 .8 ,106., -40. 

-36 .8, 106., 0. 
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PG: RIGHT ENGINE OUGHT SHE) 

6,P 

-36 .8 ,-106. ,65. 

-36 .8, -106., 0. 

-36 .8,-106. ,-40. 

-21 .8,-106., -40. 

-21 .8 ,-106., 85. 

-36 .8 ,—106 . ,85 . 

PGl LEFT ENGINE (FRCNT) 

4,P 

—21 .8 ,68. ,—40 « 

-36 .8 ,68., -40. 

-36 .8 ,106., -40. 

-21 .8 ,106., -40. 

PG: RIGHT ENGINE (FRCNT) 

4,P 

-21 .8 ,-68., -40. 

-21 .8,-106. ,-40. 

-36 .8 ,-106 .,-40 • 

-36.8 ,-68., -40. 

PG: LEFT WING (OUTER ART) 

4,P 

-36 .8, 106., 0. 

-26 .8, 213., 0. 

-26 .8 ,213., 50. 

-36 .8 ,106., 65. 

PGtRIGHT WING (OUHR ART) 

4,P 

-36 .8 ,-106., 0. 

-36 .8 ,-106. ,65. 

-26 .8 ,-213., 50. 

—26 .8 ,-213 .,0. 


RSI LEFT FUEL TANK PLATE *1 

4,P 

-26 .8 ,213., 0. 

-14.8 ,219., -32. 

-14.8 ,219., 82. 

-26 .8 ,213., 50. 

PG: RIGHT RIEL TANK PLATE #1 
4,P 

- 26 . 8 ,- 213 . , 0 . 

-26 .8,-213., 50. 

-14.8 ,-21 9., 82. 

-14.8 ,-219. ,-32. 

PG: LEFT FUEL TANK PLATE *2 
4,P 

-14.8 ,21 9., -3 2. 

—8.8 ,235. ,0 • 

-8.8 ,235., 50. 

-14.8 ,219. ,82 . 

PG: RIGHT FUEL TANK PLATE #2 

4,P 

-14.8 ,-219., -32. 

-14.8 ,-219., 82. 

-8.8 ,-23 5., 50. 

-8.8,-235.,0. 

H): ROLL ELANE 
0 . ,0 . ,90 . 

0,360,1 

F,4200. 

PP: 

T 

1,2 ^,3 

EX: 


The calculated roll plane radiation pattern shown in Figure 51 is 
found to be in good agreement with the measured pattern. 

Experimental results were obtained from NASA (Hampton, Virginia), 
using the 1/7 scale model at a range of 50 feet. 
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(b) FRONT VIEW 


Figure 49. Cessna 402B. 
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Figure 51. Roll plane pattern for Cessna 402B model. 
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Example 7: Simulation of Beechcraft Baron 


Consider a X/4 monopole mounted forward of the cockpit of a 
Beechcraft Baron aircraft as shown in Figure 52. The computer model 
based on the input data is shown in Figure 53. The input data is as 
follows: 


UN: INCHES 
3 

FQ: FIVE GIGAHERTZ 
1*5. ,1. 

FG: BEECHCRAFT BARCN 

15.5,23.2*206.5,117.5 

F 

0 ** 0 . * 0 . 

SG: MONOFQLE 

0. ,-77. 

1 

0 ** 0 . 

.414*. 828*0.*. 25*3 

1. *0. 

FG: LEFT WING 
4,T 

-10., 23 .5,-20. 

10.. 227., -10. 

10.. 227. .27. 

-10. ,23 .5*60. 

FG* LEFT ENGINE (RIGHT SIDE) 

6,F 

*-5.6,50.5,-62. 

8.4,50.5,-62. 

8.4.50.5.20. 

-8.6,50.5,20. 

-8.6,50.5,-8.5 
-5.6,50 .5,-8.5 

FG: LEFT ENGINE (TOP) 

4,F 

8.4.50.5, -62. 

11.4.84 .5, -62. 

11.4.84 .5.20. 

8.4.50.5.20. 


FG: LEFT ENGINE (FRCNT) 

4,r 

-2.6,84.5,-62. 

11*4,84 .5,-62. 

8.4,50.5,-62. 

—5.6,50.5,-62. 

FG: RIOTT WING 
4,T 

-10. ,-23.5,60. 

10., -227., 27. 

10. , -227., -10. 

-10., -23. 5, -20. 

FG: RIGHT ENGINE (LEFT SIDE) 
6,F 

-5. 6,-50. 5,-62. 

—5.6,— 50.5,— 8.5 
-8. 6,-50. 5, -8. 5 
-8.6,-50.5,20. 

8.4, -50.5,20. 

8. 4, -50. 5, -62. 

FG: RIQfT ENGINE (TOP) 

4,F 

8. 4, -50. 5,-62. 

8.4, -50.5,20. 

11.4, -84.5,20. 

11. 4, -84. 5 ,-62. 

FG: RIGHT ENGINE (RIGHT SIDE) 
6,F 

-5.6,-84.5,20. 

-5. 6,-84 .5,-7. 

-2.6,-84.5,-7. 

-2.6,-84.5,-62. 

11.4, -84.5,-62. 

11.4, -84.5,20. 
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PG: LEFT ENGINE (LEFT SIDE) 
6,F 

-5.6,84*5/20. 

11.4,84.5,20. 

11.4,84.5,-62. 

-2.6,84.5,-62. 

-2. 6, 84. 5, -7. 

-5. 6, 84. 5, -7. 


PG: RIGHT ENGINE (FRONT) 

4,F 

-2.6,-84 .5,-62. 
-5.6,-50.5,-62. 

8. 4, -50. 5, -62. 

11.4, -84.5,-62. 

H>: BOLL PLANE 
0. ,0. ,80 . 

0,360,1 

F,4200. 

PP: 

T 

1, 2.5,3 
EXs 


The conical roll plane radiation pattern is shown in Figure 54. 
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(b) FRONT VIEW 


Figure 52. Beechcraft Baron with Antenna in forward location. 
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Figure 54. Roll conical pattern (9 p =80°) for Beechcraft Baron model 
shown in Figure 53. 
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Next, let us consider the effect of the rotating propellers in the 
front of the engines. It is necessary to check the scattering due to 
the rotation of the propellers because they are close to the antenna. 
Four different positions (i.e., 0°, 45°, 90°, 135°) of the stationary 
propellers are chosen to simulate the rotating motion of the propellers 
as shown in Figure 55. For simplicity, only the left propeller is 
considered here. The input data for the four different propellers are 
as follows: 


pkupelluk (tup) At 0 
3 F 

2^9, -67. 5, -66. 

42.9, -70., -65. 

42.9, -65., -67. 

PG: PROPELLOR (BOTTOM) AT 0 
3 F 

2^9, -67. 5, -66. 

-37. 1,-70., -67. 

-37. 1,-65., -65. 


PG: PROPELLOR (TOP) AT 45 
3,F 

2. 9, -67. 5, -66. 

29.414^-97.554,-65. 

32.954,-94.014,-67. 

PG: PROPELLOR (BOTTOM) AT 45° 
3,F 

2. 9, -67. 5, -66. 

-27.154,-40.986,-67. 

-24.384,-37.446,-65. 


BG: PROPELLOR (TOP) AT 90* 
3,F 

2.9, -67. 5,-66. 
0.4,-107.5,-65. 

5.4, -107.5,-67. 

PG: PROPELLOR (BOTTOM) AT 90 
3,F 

2.9 , —67 .5,-66 . 

0. 4,-27. 5, -65. 

5. 4, -27 .5,-67. 


PG: PROPELLOR (TOP) AT 135° 
3,F 

2.9, -67. 5, -66. 

32.954^-40.986,-65. 

29.414,-37.446,-67. 

PG: PROPELLOR (BOTTOM) AT 135° 
3,F 

2. 9, -67. 5,-66. 

-23.614,-97.554,-67. 

-27.154,-94.014,-65. 
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Figure 55. Beechcraft Baron model with rotating propellers on one 
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The calculated conical roll plane radiation patterns for the four 
different propeller positions are shown in Figures 56(a) through (d). 

The roll plane pattern in Figure 57 is a combination of the four 
previous patterns. The width of the pattern line indicates the 
variation of the radiation pattern due to the rotation of the 
propellers. 
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Figure 56(a). 


Roll conical pattern for the Beechcraft Baron with 
propeller at 0° (vertical). 
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Figure 56(b). Roll conical pattern for the Beechcraft Baron with 
propeller at 45°. 
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Figure 56(c). Roll conical pattern for Beechcraft Baron with propeller 
at 90° (horizontal). 
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Figure 56(d). Roll conical pattern for the Beechcraft Baron with 
propeller at 135°. 
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Figure 57. Variation in the roll conical pattern for the Beechcraft 
Baron due to the rotation of the propellers. 
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Example 8: Simulation of the Cessna 150 Aircraft 


In this example, a monopole antenna is mounted forward on the wings 
of a Cessna 150 aircraft. A different approach was taken in the 
modeling of the Cessna 150 aircraft due to the position of the antenna 
and the shape of the wings. The wings are simulated by the composite 
ellipsoid and two flat plates are attached to simulate the nose and 
fuselage of the aircraft. The line drawing of the Cessna 150 is shown 
in Figure 58, and the computer model based on the input data is shown in 
Figure 59. The input data is as follows: 


UN: INCHES 
3 


cn. 
• y* 


1.5. ,1. 

FG: CESSNA 150 WING 
3.5,25., 250. ,250. 

F 

0 ., 0 ., 0 . 

FC: SQUARE OFF WING TIPS 


T,T 

196.25,-196.25 

SG: MONOPOLE MOUNTED ON WING 

66 . , 0 . 

1 


o.,o. 

.414, .828,0. ,.25,3 

1 ., 0 . 

PG: NOSE PLATE 
4,T 

-17.35,32.53,-17.52 

-27.35,82.4,-15.18 

-27.35,82.4,15.18 

-17.35,32.53,17.52 


PG: FUSELAGE PLATE 
4,T 

-18. 67 ,-25 ,2,1 4. 93 

10., -152.,3.17 

10. , -152. ,-3.17 
-18.67,-25.2,-14.93 

PD: ELEVATION PLANE AS REFERED TO THIS MODEL 
0. , 0 . ,90 . 

0,360,1 

F,4200 

PP: PEN PLOT 
T 

1,1.925,3 

EX: 


The resulting elevation plane pattern for this model is shown in 
Figure 60. Although the magnitude of the ripple is not quite perfect, 
it is of the correct spatial frequency, and the general shape of the 
pattern is good. 
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Figure 58. Cessna 150 




(a) Top view. 
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(b) Side view. 


Figure 59. Cessna 150 model. Dashed lines are not part of the computer 
simulation. 
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Figure 60. Elevation plane pattern for a Cessna 150 aircraft. 
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Example 9: Simulation of F-16 Fighter Aircraft 


Consider a TACAN antenna mounted on the top of a F-16 fighter 
aircraft as shown in Figure 61 and operated at a frequency of 0.96 GHz. 
A composite ellipsoid and a total of 12 plates are used to simulate the 
structure of the aircraft. The computer model is illustrated in Figure 
62. Note that the radome of the F-16 is simulated as a truncated 
fuselage. The measured data was obtained by General Dynamics using a 
quarter scale model of the F-16. The input data is as follows: 


UN: IN INCHES 
3 

FG:F16A FUSELAGE GEOMETRY AT STATION 250 
*>i « ox ^nn ?*n 

4. I • o' f | -TW W • I V m 

F 

o.,o.,o. 

FC: 

T.T 

300.. -185. 

FQ: FREQUENCY 

1.0. 96.1 . 

SG: SOURCE GEOMETRY 
0., 13.25 
1 

o.,o. 

0.,0.,0., 3.07 58,3 

1 .. 0 . 

PGjCURVATURE SIMULATED PLATE #1 ON P0S. SIDE 
6,T 

8.2046,22.4421,-61 . 

2.1418,36.5,-61 . 

-4.0866,50.942,-8.6 
-5.4054, 54. ,8.743 
-5. 4054, 54., 158 .95 
8.2046,22.4421,158.95 

PG: CURVATURE SIMULATED PLATE #2 ON POS. SIDE 
5,T 

8.2046,22.4421 ,158.95 
-5. 4054, 54., 158. 95 
-6. 2805, 54., 209. 084 
-7 .6944,54. ,290. 084 
5.9156,22.4421,290.084 


PG: CURVATURE SIMULATE 

PLATE #3 ON POS. SIDE 


3 T 

o!5,19.2,-150. 

2.1418,36.5,-61. 

8.2046,22.4421,-61. 
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PGjWING on pos. side 
4,F 

-5. 4054,54., 8. 743 
-5.4054, 180. ,114.47 
-5 .4054, 180., 158. 95 
-5. 4054, 54., 158 .95 

PG: HORIZONTAL STABILIZER ON POS. SIDE 
4,F 

-5.4054,54. ,21 9.7958 

-5.4054,109.101 ,266.031 

-5.4054,109.101 ,290.084 

-5. 4054,54., 290. 084 

PG: VERTICAL STABILIZER ON NEG. SIDE 

4»T 

20 ., 0 ., 1 60 . 

120. ,0., 261. 

120.. -3 .4, 298. 

20.. -6 .8,234. 

PG: VERTICAL STABILIZER ON POS. SIDE 

4, T 

20. . 6. 8.234. 

120. . 3.4. 298. 

120. . 0. .261 . 

20 . , 0 . , 1 60 . 

PG: CURVATURE SIMULATED PLATE #1 ON NEG. SIDE 
6,T 

8.2046, -22.4421 ,158.95 
-5 .4054, -54., 158 .95 
-5. 4054, -54., 8. 743 
-4.0866,-50.942,-8.6 
2.1418,-36.5,-61 . 

8.2046, -22.4421,-61. 

PG: CURVATURE SIMULATED PLATE #2 ON NEG. SIDE 

5, T 

5.9156,-22.4421 ,290.084 
-7 .6944, -54., 290. 084 
-6 .2805 ,-54., 209. 084 
-5 .4054 ,-54., 158 .95 

8.2046, -22.4421,158.95 

PG: CURVATURE SIMULATED PLATE #3 ON NEG. SIDE 

3, T 

8.2046, -22.4421,-61. 

2.1418,-36.5,-61. 

0.5,-19.2,-150. 

PG:WING ON f€G. SIDE 

4, F 

-5 .4054, -54., 158 .95 
-5. 4054, -180., 158 .95 
-5. 4054 ,-180., 114.47 
-5 .4054,-54., 8 .743 
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PG: HORIZONTAL STABILIZER ON NEG. SIDE 
4,F 

-5 .4054, -54., 290. 084 
-5.4054,-109.101 ,290.084 
-5.4054,-109.101 ,266.031 
-5.4054,-54., 21 9.7958 
PP: POLAR PLOT IN DB 
T 

1,2.81,3 

PD: AZIMUTH PLANE CUT 

90.,0.,10. 

0,360., 1 
T, 50000. 

EX: EXECUTE 

PD: AZIMUTH PLANE CUT 

90 . , 0 . ,20 . 

0,360., 1 
T,50000. 

EX: EXECUTE 

PD: AZIMUTH PLANE CUT 

90 .. 0..30. 

0,360. ,1 
T,50000. 

EX: EXECUTE 
I 

I 

1 

PD: AZIMUTH PLANE CUT 

90. . 0.. 115. 

0,360. ,1 

T, 50000. ' 

EX: EXECUTE 

PD: AZIMUTH PLANE CUT 

90. . 0.. 120. 

0,360. ,1 

T, 50000. 

EX: EXECUTE 


PD:ELEVATI0N PLANE CUT 
90. ,90. ,90. 

0,360. ,1 
T,50000. 

EX: EXECUTE 
PI: 

9 

PD: ROLL PLANE CUT 
0. ,0 . ,90 . 

0,1 80., 1 
T,50000. 

TO: 

F,F,F 

T,T 

T,T,T,T,T,T,T,T 

1.9.1 

2 . 6.1 

1.5.1 

1.3.2 
1,4,1 
1,4,1 
1,4,1 

1.4.1 

2.3.1 

1.3.2 
DD: 

5 

1,6,1, 3 

1. 6.4.1 

1.6 ,4 ,4 

4.1 ,1,6 

4.4. 1.6 
RD: 

1 

1,4 

EX: EXECUTE 


To show the complete volumetric radiation patterns, the various 
azimuthal conical patterns are calculated as shown in Figures 63 through 
82. In each case, both the principal and cross polarizations are 
considered. The elevation plane and roll plane patterns are also shown 
in Figure 83 and 84, respectively. All the above calculated results 
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compare favorably with the measurements. It is noted that since the 
cockpit section simulation is not complete in our model, one cannot 
expect good agreement between the calculated and measured results in the 
nose region since the cockpit is part of the radiation path. In 
addition, the ripple above the aircraft in the elevation pattern are 
most likely created by the cockpit which is not simulated in this 
model . 
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(c) TOP VIEW 


Figure 62. Computer simulated model of a F-16 fighter aircraft. 
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Figure 63. Azimuthal conical pattern (0p=io°) of a A/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 64. Azimuthal conical pattern (Qp=20°) of a x/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 65. Azimuthal conical pattern (9p=30°) of a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 66. Azimuthal conical pattern (Qp=40°) of a x/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 67. Azimuthal conical pattern (9p=45°) 0 f a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 68. Azimuthal conical pattern (9p=50°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 69. Azimuthal conical pattern (0p=55°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 72. Azimuthal conical pattern (9n=70°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 71. Azimuthal conical pattern (9p=65°) of a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 70. Azimuthal conical pattern (9p=60°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 73. Azimuthal conical pattern (9p=75°) of a \/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 74. Azimuthal conical pattern (9 p =80°) of a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 75. Azimuthal conical pattern (9n=85°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 76. Azimuthal conical pattern (Qp=90°) of a x/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 77. Azimuthal conical pattern (9p=g5°) 0 f a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 78. Azimuthal conical pattern (9p=l00°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 79. Azimuthal conical pattern (9p=l05°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 80. Azimuthal conical pattern (9p=no°) of a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 81. Azimuthal conical pattern (9p=ll5°) of a X/4 monopole 
mounted on top of a F-16 fighter aircraft. 
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Figure 82. Azimuthal conical pattern (9p=l20°) of a X/4 monopole 
mounted on top of a F-16 fignter aircraft. 
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Figure 83. Elevation plane pattern of a X/4 monopole mounted on top of 
a F-16 fighter. 
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Figure 84. Roll plane pattern of a x/4 monopole mounted on top of a 
F-16 fighter aircraft. 
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Example 10: Simulation of F-4 Fighter Aircraft 


Consider a x/4 monopole mounted on the bottom fuselage of an F-4 
aircraft, which is loaded with armament, and operated at a frequency of 
.375 MHz. The measured data was obtained at the RADC Newport antenna 
range. The line drawings and the computer model of the F-4 aircraft are 
illustrated in Figure 85 and 86, respectively. Note that since the 
antenna is mounted on the belly of an aircraft, the coordinates are 
defined so as to associate with the bottom part of the aircraft. 
Consequently, the geometry of the F-4 in our computer model, as well as 
the pattern coordinate systems, are turned upside down. In fact, for 
the 9p=75° pattern computed here, the actual pattern angle from the 
vertical is 180° - 75° or 105°. The input data is as follows: 


UN: INCHES 
3 

PC: FREQUENCY 375 WZ 

1. . 375.1. 

PG: F-4 FUSELAGE 

5.. 20..300..250. 

F 

0 . , 0 . , 0 . 

SG: MONOBOLE 

0 .,- 200 . 

1 

0. fO. 

.414,. 828,0. ,7.87,3 

1.. 0. 

PG: LEFT WING 
6,T 

-2. ,18., -133. 

-2. ,50., -133. 

-2. ,50. ,-70. 

-2., 230., 119. 

-2., 230., 167. 

-2. ,18., 136. 


PG: RIGHT WING 
6,T 

-2., -18. ,136. 

-2. ,-230., 167. 

-2. ,-230., 119. 

—2. ,—50 . ,—70. 

-2., -50., -133. 

-2., -18., -133. 

FGiLEFT ENGINE INTAKE 

4,T 

-2., 18., -133. 

-2., 50., -133. 

-4., 50., -133. 

-4., 18., -133. 

PG: RIGHT ENGINE INTAKE 
4,T 

-4., -18., -133. 

-4. ,-50., -133. 

-2. ,-50., -133. 

-2., -18., -133. 
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PG: FUSELAGE BLOCKAGE 

«/r 

4. . 9. ,-50. 

15.. 13. ,-50. 

25. ,6.,— 50". 

25., -6.,-50. 

15 • 13 . , — 50 . 

4., -9.,-50. 

FGsLEET engine intake 
4,F 

-4. ,50., -133. 

-52., 50., -133. 

-52. ,18. ,-133. 

-4. ,18., -133. 

PG: RIGHT ENGINE INTAKE 
4,F 

-4.,-18.,-133. 

-52., -18., -133. 

-52. ,-50. ,-133. 

-4. ,-50., -133. 

PG: LEFT MISSILE PYLCN 
4,F 

20. . 72. ,-45. 

0.,72.,-45. 

0.,72.,-117. 

20.. 72., -117. 


PG: RIGHT MISSILE PH/M 

4,F 

20., — 72.,— 117 • 

0., -72., -117. 

0.,— 72.,— 45. 

20., -72., -45. 

PGxLEFT fuel tank 
4,F 

36. . 127. .148, 

0.,127.,148. 

0.,127.,-77. 

36. . 127. ,-77. 

PG: RIGHT FUEL TANK 
4,F 

36., -127., -77. 

0., -127;, -77. 

0., -127., 148. 

36., -127., 148. 

PP: POLAR PLOT 
T 

1,2.9, 3 

PD:AZM7IH CONICAL PATTERN 

90. . 0. .75, 

0,360,1 

F,5000. 

Bit 


The azimuthal conical pattern is compared with measured data in 
Figure 87. Although there exists some discrepancy, the general shape of 
the two patterns are in good agreement. 
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(a) SIDE VIEW 



Figure 85. F-4 (Phantom) fighter aircraft. 
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Figure 86. Computer simulated model of a F-4 Phantom fighter aircraft. 
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Figure 87. Azimuthal conical pattern (00*105°) of a A/4 monopole 
mounted on the belly of a F-4 fighter aircraft. 
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Example 11: Simulation of an A-10 Aircraft 

Consider four monopoles mounted on the belly of an A-10 aircraft as 
shown in Figure 88. Each monopole is spaced a half wavelength apart. 

The mutual coupling between the monopoles is significant and cannot be 
neglected in the pattern calculations. The excitation Including the 
coupling effect on the radiators is obtained using the thin-wire moment 
method [4]. The computer model based on the input data is shown in 
Figure 89. The input data is as follows: 


UN: INCHES 

3 

FQ: 17.576 (HZ 

1,17.576,1 

PG: A-10 FUSELAGE 

0.117,0.425,8.05,2.62 

F 

Sg!°FCWOR 3LE ARRAY WITH (DUELING EFFECT INCLUDED IN CURRENT VALUE 
O.i-1.29 

4 

.336,180. 

0. ,0. ,0. ,0.168,3 
.272,14. 

.336,0. 

0.,0.,0. ,0.168,3 
.272,14. 

.168,270. 

0.,0.,0., 0.168,3 

1.0,-3.0 

.504,270. 

0 . ,0. ,0 . ,0 .168,3 
.272,14. 

PG: LEFT WING INNER PLATE 
4,T 

0.05,. 425,. 07 
0.05, 1.86, .07 
0.05,1.86,2. ' 

0.05, .425,2. 

PG: RIGHT WING INNER PLATE 
4,T 

0.05, -.425, 2. 

0.05,-1.86,2. 

0.05,-1.86,. 07 
0.05,-.425,.07 
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PG: LETT WING OUIBt IUB 
4,F 

0.05, 1.86, .07 
-.45, 5. 6, .49 
-0.45,5.6,1.77 
0.05,1.86,2. 

FG: RIGHT WING COT® HATE 
4,F 

0.05,-1.86,2. 
-.45,-5.6,1.77 
-.45,-5. 6, .49 
0.05 ,-1.86, .07 
FG: LEFT FUELTANK 
4,F 

0.05, 1.58, -.6 

•30,1.58,-.6 

.30,1.58,2. 

0.05,1.58,2. 

FG:LEET FUEL TANK 
4,F 

•30, 1.58, -.6 
.59,1.86 ,-.6 
.59,1.86,2. 

.30,1.58,2. 

FG: RIGHT FUEL TANK 
4,F 

.30,-1.58, 2. 

.59,-1.86,2. 

.59,-1.86 ,—.6 
.30,-1.58,-.6 
FG: RIGHT FUEL TANK 
4,F 

0.05,-1.58,2. 

.30,-1.58,2. 

.30,— 1.58,— .6 
0.05,-1.58,-.6 


FG: KLON 1 
4,F 

•4,0.53,1.61 
.05,0.53,1.61 
•05,0.53, .56 
•4,0.53, .33 
FG: KLCN 2 
4,F 

.4, -.53,. 33 
•05,-.53,.56 
.05 ,-.53 ,1.61 
•4,-.53,1.61 
FG: fxlcn 3 
4,F 

•4,1.05,1.61 
.05,1.05,1.61 
•05,1.05, .56 
•4,1.05, .3 
FG: FXLCN 4 
4,F 

.4.-1. 05. .3 
.05,-1.05, .56 
.05,-1.05,1.61 
.4,-1.05,1.61 

PD: AZDtTEH CONICAL PATTERN 
90.,0.,75. 

0,360,1 

T,6000. 

FP: 

T 

1,2.5, 3 
EK: 


The azimuthal conical patern (9p*105°) is compared with the 
measured data obtained at the RADC Newport site in Figure 90. The 
result shows good engineering agreement. 
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(c) TOP VIEW 


Figure 88. A-10 aircraft. 
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Figure 90. Azimuthal conical pattern (9p=l05°) of four monopoles 
mounted on the belly of a A-IO aircraft. 
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Example 12: Simulation of C-141 Aircraft 

Consider a monopole mounted on the top of a C-141 aircraft. The 
line drawings and the computer model of the C-141 are shown in Figures 
91 and 92, respectively. The input data is as follows: 


UN: PEET 
2 

FG: C141 FUSELAGE GEOMETRY 
7. 37,8. 37, 90. ,46. 05 
F 

0. ,0 • ,0 • 

FQ: 2.52 GHZ 

1. ,2.52,1. 

SG: SOURCE GEOMETRY 

0 . , 0 . 

1 
X 

0 . r 0 • 

0.2, 0.5,0. ,0.09, 3 

1.. 0. 

PG: WING ON POSITIVE SIDE 
5 T 

6^0, 7. 37 ,-7.4 

6.0. 78.6.27.04 

6.0. 78.6.36.86 

6.0. 30.7.24.6 

6.0. 7.37.23.1 
PG: WING ON NEGATIVE SIDE 
5 T 

6 !o, -7. 37, 23.1 

6.0, -30.7,24.6 

6.0, -78.6,36.86 

6.0, -78.6,27.04 

6. 0, -7. 37 ,-7.4 
PG: VERTICAL STABILIZER POSITIVE SIDE 
4 ,T 

7.5,1.6465,88.11 

24. 58.1. . 92. 81 

24.58.0. ,77 .45 

7. 5.0. .62. 82 

PG: VERTICAL STABILIZER NEGATIVE SIDE 
4 ,T 

7. 5.0. .62. 82 

24. 58.0. .77 .45 
24. 58.-1. ,92. 81 
7.5,-1.6465,88.11 


FG: T-TAIL POSITIVE SIDE 

4,F 

24. 58.1. . 92.81 
24.58.25.3,98.22 
24.58.25.3,92.05 

24. 58.0. .77 .45 

PG: T-TAIL NEGATIVE SIDE 
4,F 

24. 58,0. ,77 .45 

24.58. -25.3,92.05 

24.58, -25.3,98.22 
24. 58-1. ,92. 81 

PP: POLAR PLOT IN DB 
T 

1,1.42,3 

PD: AZIMUTH PLANE PATTERN 
90 . ,0 . ,90 . 

0,360,1 

F,1000. 

EX: 
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Figure 91. C-141 aircraft. 


2i 





(o) SIDE VIEW 
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Figure 92. Computer simulated model of a C-141 aircraft. 
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Various azimuthal conical patterns (10°, 20°, 30°, 40°, 50°, 60°, 
70°, 80°, 90°, 100°, see Figure 10(c)) and the elevation plane pattern 
are computed and compared with measured results as shown in Figures 93 
and 94, respectively. The experimental work was performed at General 
Dynamics (San Diego, California) using a 1/10 scale model of a C-141 
aircraft. The calculated results compare very favorably with the scale 
model measurements. 
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Figure 93. Azimuthal conical patterns of a X/4 monopole mounted on a 
C-141 aircraft. 
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Figure 93. (Continued). 
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Figure 93. (Continued). 
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Figure 93. (Continued). 
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Figure 93. (Continued). 
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Figure 94. Elevation plane pattern of a x/4 monopole mounted on top of 
a C-141 aircraft. 
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Example 13: Simulation of Missiles 


Consider an axial slot antenna mounted between two large ram jets 
on a missile. The front view of the missile and the computer model are 
shown in Figure 95. The fuselage is simulated as a composite ellipsoid 
and only two side walls of the ram jets are simulated with two flat 
plates. A fuselage blockage plate is also added in the computer model 
and is shown as a dash line in the corresponding figure. This feature 
is added in the code to prevent the rays from passing through the 
fuselage and contributing to the field calculations. The roll and 
elevation plane patterns are compared with the measured results in 
Figure 96 and 97, respectively. Again, the calculated and measured 
results are in good agreement. The input data is as follows: 


UN: IN INCHES 

3 

F6: FUSELAGE GEOMETRY FOR THE ANTENNA MOUNTED BETWEEN TWO RAM JETS CASE 

6 . 75 . 6 . 75 . 100 .. 100 . 

F 

0 «, 0 «, 0 « 

FC: 

T,T 

21 ., - 21 . 

FQ: FREQUENCY 

1 , 12 . , 1 . 

SG: SOURCE 

1 . 32 . 0 . 

1 

0 . , 0 . 

0 . 024 , 0 . 153 , 0 . , 0 . 075,1 

1 .. 0 . 

FB: FUSELAGE BLOCKAGE 
1 

4 

6 . 1 , 8 . 3 , - 20 . 

6 . 1 , 8 . 3 , 20 . 

6 . 1 , - 8 . 3 , 20 . 

6 . 1 , - 8 .3 ,- 20 . 

/ 
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PG: RIGHT AIR DUCT #1 
4,T 

6.147,2.788,-20. 

7.826,4.467,-20. 

7.826.4.467.10. 

6.147.2.788.10. 

PG: RIGHT AIR DUCT #2 
4#T 

6.147, -2.788,10. 

7.826, -4.467,10. 

7.826, -4.467,-20. 

6.147, -2.788,-20. 

PG: LEFT AIR DUCT #1 
4.T 

6.147.2.768.10. 

7.826.4.467.10. 

7.826.4.467.20. 

6.147.2.788.20. 

PG: LEFT AIR DUCT #2 
4»T 

6.147, -2.788,20. 

7.826, -4.467,20. 

7.826, -4.467,10. 

6.147, -2.788,10. 

PO: ROLL PLANE PATTERN 
0., 1.32,90. 

0,360,1 

F,600. 

PP: PEN PLOT 
T 

1 .3.. 3 

EX: EXECUTE 

H): ELEVATION PLANE PATTERN 

90. . 91 .32.90.2 
0,360,1 
F,600. 

EX: EXECUTE 
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(a) MISSILE MODEL 



(b) COMPUTER SIMULATED MODEL 


Figure 95. Missile model for an axial slot mounted between two ram 
jets. 
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Figure 96. Roll plane pattern for an axial slot mounted between two 
ram jets. 
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Figure 97. Elevation plane pattern for an axial slot mounted between 
two ram jets. 
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In this case, an axial slot antenna Is mounted on the ram jet instead of 
the fuselage as shown in Figure 98. The ram jet is simulated as a 
composite ellipsoid and the other structures are simulated by multiple 
flat plates. The roll and elevation plane patterns are compared with 
measured data In Figure 99 and 100, respectively. The input data Is as 
follows: 


UN: IN INCHES 
3 

F6: FUSELAGE GEOMETRY FOR THE ANTENNA 

2.375.2.375.100. . 100. 

F 

o.,o.,o. 

FC: 

T,T 

16. , -16. 

FQ: FREQUENCY 

1 , 12 . , 1 • 

SG: SOURCE GEOMETRY 

o.,o. 

1 


0 . , 0 . 

0.024,0. 153,0. ,0.075,1 

1 ., 0 . 

PG: LEFT SIDE WALL 
4,T 

-0.3,-2.3,15. 

-1.9541 ,-2.375,15. 
-1.9541 ,-2.375,-15. 
-0.3,-2.3,-15. 

PG: LEFT SIDE SLOPE 
4,F 

-1.9541 ,-2.375,15. 
-6.75,-7.1709,15. 
-6.75,-7.1709,-15. 
-1.9541 ,-2.375,-15. 

PG: LEFT AIR DUCT 
4,F 

-6.75,-7.1709,15. 

-6.75,-9.5,15. 

-6. 75, -9. 5, -15. 
-6.75,-7.1709,-15. 

PG: RIGHT SIDE WALL 
4,T 

-0.3 ,2. 3, -15. 

-1.9541 ,2.375,-15. 

-1 .9541 ,2.375,15. 
-0.3,2.3,15. 
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MOUNTED ON THE RAM JET CASE 


PG: RIGHT SIDE SLOPE 
4,F 

-1.9541 ,2.375,-15. 

-6.75,7.1709,-15. 

-6.75,7.1709,15. 

-1.9541 ,2.375,15. 

PG: RIGHT AIR DUCT 
4,F 

-6.75,7.1709,-15. 

-6.75,9.5,-15. 

-6.75,9.5,15. 

-6.75,7.1709,15. 

PD: ROLL PLANE PATTERN 
0. ,0. ,90 . 

0,360,1 

F,600. 

PP: PEN PLOT 
T ■ 

1 *3. ,3 
EX: EXECUTE 

H>: ELEVATION PLANE PATTERN 
90. ,93.76,90. 

0,360,1 

F,600. 

EX: EXECUTE 





(b) COMPUTER SIMULATED MODEL 


Figure 98. Missile model for an axial slot mounted on a ram jet. 
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Figure 99. Roll plane pattern for an axial slot mounted on a ram jet. 
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Figure 100. Elevation plane pattern for an axial slot mounted on a ram 
jet. 
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Example 14: Simulation of Space Shuttle 


Consider an S-band Quad antenna (crossed-slot antenna with a 90° 
phase difference between the two slots) mounted on the top of a Space 
Shuttle Orbiter as shown in Figure 101. The computer model of the Space 
Shuttle Orbiter based on the input data is illustrated in Figure 102. 

The three principal plane patterns and a 45° roll conical pattern are 
computed as shown in Figures 103 to 106. The input data is as follows: 


UN: IN INCHES 
3 

FG: SPACE SHUTTLE FOR THE TOP MOUNTED ANTENNA KITH CLOSED PAYLOAD DOORS CASE 

75. . 104. .1500.. 130. 

F 

0 . , 0 . , 0 . 

FC: 

T#F 

944.3.0. 

FQjFREQUENCY (2 GHZ) 

1 . 2 .. 1 . 

SG: CROSSED SL0T(.39 *.78 WAVELENGTH) 

62.5,-10. 

2 

0 ., 0 . 

2.3.4.6.0. . 1.476.1 

1 .. 0 . 

0 ., 0 . 

2.3.4.6.90. . 1.476.1 

1.. 90. 

PGjSIDE WALL SIMULATED PLATE #1 ON POS. Y-SIDE 

4.T 

0., 101. 1928,-30. 

-120. ,102.3928,-30. 

-109.35,110.58,100. 

-5., 109. 54, 100. 

PG.-SIDE WALL SIMULATED PLATE #2 ON POS. Y-SIDE 
4,T 

-5., 109. 54, 100. 

-109.35,110.58,100. 

-78.77,117.94,381. 

-5. ,117.2,381. 
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PGsSIDE HALL SIMULATED PLATE #3 ON POS. Y-SIDE 
4,T 

-5.,1 17.2,361. 

-78.77,1 17 .94*381 • 

-120. ,118.84,458. 13 
-5.,1 17.69,458.13 

PGjSIDE HALL SIMULATED PLATE #4 ON POS. Y-SIDE 
4»T 

-5. ,117.69,458.13 
-120. ,118.84,458.13 
-147..1 06.07 ,944.3 
-5., 104 .65, 944 .3 

PGrWING SIMULATED PLATE #1 ON POS. Y-SIDE 

3. F 

-120., 102.3928,-30. 

-120. ,125. 192, 100. 

-109.35,110.58,100. 

PGjHING SIMULATED PLATE #2 ON POS. Y-SIDE 

4. F 

-109.35,110.58,100. 

- 1 on 1 ok 109 inn 

"i4U«f • 4./ • • #*> f I wv# 

-120. ,188. ,458. 13 
-78.77,117.94,381 . 

pGsWING SIMULATED PLATE #3 ON POS. Y-SIDE 

3, F 

-78.77,117.94,381. 

-120. ,188. ,458. 13 
-120. ,118.84,458.13 

PG.-WING SIMULATED PLATE #4 ON POS. Y-SIDE 

5. F 

-120. .118.84.458.13 
-120. ,188. ,458. 13 
-135.083,468.34,729.72 
-142.739,468.34,867.57 
-147., 106 .07 ,944.3 

PG:VERTICAL STABILIZER ON POS. Y-SIDE 

4, T 

83. . 14. 1.884.5 
390.7,5.8,1098.7 

390 .7.0. . 1033 .6 

83. . 0. .725.9 

PG: VERT I CAL STABILIZER ON NEG. Y-SIDE 
4,T 

83. . 0.. 725.9 

390 .7.0. , 1033.6 
390.7,-5.8,1098.7 

83. , -14.1 ,884 .5 
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PGsNOSE SIMULATED PLATE ON POS. Y-SIDE 
4,F 

-65. ,50. ,-322. 

-105., 50. 4, -322. 

-120., 102.3928,-30. 

0., 101. 1928,-30. 

PGsNOSE SIMULATED PLATE #1 ON TOP 
4,T 

-22.26,0. ,-130. 

-65., 0., -322. 

-65., 50., -322. 

0.,101 .1928,-30. 

PGsNOSE SIMULATED PLATE #2 ON TOP 

4, T 

0., -101 .1928,-30. 

-65. ,-50. ,-322. 

-65., 0., -322. 

-22.26,0. ,-130. 

PGsNOSE SIMULATED PLATE ON FRONT 

5, F 

-105., 50. 4, -322. 

-65. ,50. ,-322. 

-65., 0., -322. 

-65., -50., -322. 

-105. ,-50.4,-322. 

PGsNOSE SIMULATED PLATE ON NEG. Y-SIDE 
4,F 

0.,-101 .1928,-30. 

-120., -102. 3928, -30. 

-105. ,-50 .4 ,-322. 

-65., -50., -322. 

PGsSIDE NALL SIMULATED PLATE #1 ON NEG. Y-SIDE 

4,T 

-5. ,-109.54,103. 

-109.35,-110.58,100. 

-120. ,-102.3928,-30. 

0., -101. 1928,-30. 

PGsSIDE NALL SIMULATED PLATE #2 ON NEG. Y-SIDE 

4.T 

-5. ,-117.2,381. 

-78.77,-117.94,381. 

-109.35,-110.58,100. 

-5., -109. 54, 100. 

PGsSIDE NALL SIMULATED PLATE #3 ON NEG. Y-SIDE 
4,T 

-5., -117.69,458.13 
-120. ,-118.84,458. 13 
-78.77,-117.94,381. 

-5.,-1 17. 2,381. 


222 



PGjSIDE WALL SIMULATED PLATE #4 ON NEG. Y-SIDE 
4*T 

-5., -104. 65, 944 .3 
-147., -106 .07 ,944. 3 
-120. ,-118.84 ,458. 13 

PGjWING SIMULATED PLATE #1 ON NEG. Y-SIDE 

3, F 

-109.35,-110.58,100. 

-120. ,-125. 192, 100. 

-120., -102. 3928, -30. 

PG:WING SIMULATED PLATE #2 ON NEG. Y-SIDE 

4, F 

-78.77,-117.94,381. 

-120., -188., 458. 13 
-120. ,-125.192, 100. 

-109.35,-110.58,100. 

PG:WING SIMULATED PLATE #3 ON NEG. Y-SIDE 


3.F 

-120. ,-118.84,458. 13 
-120. ,-188. ,458. 13 


no nn . 
i i $' 


.in n j ifii 

'II# 


PG:WING SIMULATED PLATE #4 ON NEG. Y-SIDE 


5,F 

-147., -106. 07, 944 .3 
-1 42 .739 ,-468 .34 ,867 .57 
-135.083 ,-468 .34,729.72 
-120. ,-188. ,458. 13 
-120. ,-118.64,458. 13 
PP: POLAR PLOT IN DB 
T 

1,2.766,3 


PDsROLL CONICAL PLANE CUT 
0 . ,0. ,45 . 

0,360,1 
T, 50000. 

EX: EXECUTE 


PD: ELEVATION PLANE CUT 

90. . 90. .90. 

0,360,1 

T, 50000. 

EX: EXECUTE 

PD: AZIMUTH PLANE CUT 

90.. 0..90. 

0,360,1 

T, 50000. 

EX: EXECUTE 

PD: ROLL PATTERN CUT 

0.,0.,90. 

0,360,1 
T, 50000. 

EX: EXECUTE 
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(a) 3-DIMENSIONAL VIEW 




(d) TOP VIEW 


Figure 102. Computer simulated model for a crossed-slot antenna mounted 
on top of a Space Shuttle Orbiter. 
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Figure 103. Elevation plane patterns for a crossed-slot antenna mounted 
on top of a Space Shuttle Orbiter. 
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Figure 104. Azimuth plane patterns for a crossed-slot antenna mounted 
on top of a Space Shuttle Orbiter. 
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(SCALE: EACH DIVISION** 40B) 
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Figure 105. Roll plane radiation patterns for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter. 
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Figure 106. Roll conical patterns (0 p =45°) for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter when the payload 
bay doors are closed . 
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When the payload doors of the Orbiter are open, six more plates are 
used to simulate the heat radiators inside the payload doors, but at the 
same time ten plates are removed from the model since these plates are 
in the shadow region of the payload doors. The computer model of the 
Orbiter with open payload doors is shown in Figure 107. 

In order to analyze the blocking and reflection effect of the heat 
radiators, the roll conical patterns (9p=45°) are computed and then 
compared with the roll conical patterns obtained when the payload doors 
were closed. The input data for this case is as follows: 


UN: IN INCHES 

3 

F6: SPACE SHUTTLE FOR THE TOP MOUNTED ANTENNA WITH OPEN PAYLOAD DOORS CASE 

75. . 104. .1500. .130. 

F 

0 . , 0 . , 0 . 

FC: 

T,F 

19. . 0. 

FQ: FREQUENCY (2 GHZ) 

1 . 2 .. 1 . 

SGjCROSSED SLOT ( .39 *.78 WAVELENGTH) 

62.5,-10. 

2 

0 ., 0 . 

2.3.4.6.0. . 1.476.1 

1 .. 0 . 

0 ., 0 . 

2.3.4.6.90. . 1.476.1 

1 .. 90. 

FB: 

2 

4 

-5., 104 .37, 19. 

-5. ,104 .37,720. 

-5., -104 .37, 720. 

-5. ,-104 .37, 19. 

4 

83.. 104 .37. 720. 

83., -104 .37, 720. 

-83. ,-104.37,720. 

-83. ,104 .37,720. 
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PGjSIDE MALL SIMULATED PLATE #1 ON POS. Y-SIDE 
5,T 

0., 101. 1928,-30. 

-120., 102.3928,-30. 

-109.35,110.58,100. 

-5. ,109 .54, 100. 

-5. ,104 .37, 19. 

PG:DOOR #1 ON POS. Y-SIDE 
4,F 

-5. ,104 .37, 19. 

-37.1,130.1,19. 

-37.1,130.1,720. 

-5., 104 .37 ,720. 

PGsDOOR #2 ON POS. Y-SIDE 
4,F 

-37.1,130.1,19. 

-37.1,187.9,19. 

-37.1,187.9,720. 

-37.1,130.1,720. 

PGsDOOR #3 ON POS. Y-SIDE 

4, F 

-37.1,187.9.19. 

-5. ,245 .8, 19. 

-5. ,245 .8,720. 

-37.1,187.9,720. 

PGsWING SIMULATED PLATE #1 ON POS. Y-SIDE 

3, F 

-120., 102. 3928 ,-30. 

-120. ,125. 192, 100. 

-109.35,110.58,100. 

PGsWING SIMULATED PLATE #4 ON POS. Y-SIDE 

5, F 

-120. ,118. 84,458. 13 
-120. ,188. ,458. 13 
-135.083,468.34,729.72 
-142.739,468.34,867.57 
-147., 106. 07 ,944 .3 

PG: VERT I CAL STABILIZER ON POS. Y-SIDE 

4, F 

83. . 14. 1.884.5 
390.7,5.8,1098.7 

390. 7.0. . 1033 .6 

83. . 0. .725.9 

PG: VERTICAL STABILIZER ON NEC. Y-SIDE 
4,F 

83. . 0. .725.9 

390. 7.0. , 1033.6 
390.7,-5.8,1098.7 

83., -14. 1,884 .5 
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PGsNOSE SIMULATED PLATE ON POS. Y-SIDE 
4,F 

-65., 50., -322. 

-105., 50 .4 ,-322. 

-120., 102.3928,-30. 

0., 101. 1928,-30. 

PGsNOSE SIMULATED PLATE #1 ON TOP 

4,T 

-22. 26,0., -130. 

-65., 0., -322. 

-65., 50., -322. 

0., 101. 1928,-30. 

PGsNOSE SIMULATED PLATE §2 ON TOP 

4, T 

0., -101. 1928,-30. 

-65., -50., -322. 

-65. ,0., -322. 

-22.26,0. ,-130. 

PGsNOSE SIMULATED PLATE ON FRONT 

5, F 

-105., 50 .4, -322. 

-65., 50., -322. 

-65., 0., -322. 

-65., -50., -322. 

-105., -50. 4, -322. 

PGsNOSE SIMULATED PLATE ON NEG. Y-SIDE 

4, F 

0., -101. 1928,-30. 

-120. ,-102.3928,-30. 

-105. ,-50.4,-322. 

-65., -50., -322. 

PGsSIDE WALL SIMULATED PLATE t\ ON NEG. Y-SIDE 

5. T 

-5. ,-104 .37, 19. 

-5. ,-109.54,100. 

-109.35,-110.58,100. 

-120., -102.3928,-30. 

0., -101. 1928,-30. 

PGsDOOR #1 ON NEG. Y-SIDE 
4.F 

-5., -104 .37, 720. 

-37.1,-130.1,720. 

-37.1,-130.1,19. 

-5. ,-104 .37, 19. 
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PGjDOOR #2 ON NE6. Y-SIDE 
4 ,F 

-37.1,-130.1,720. 

-37.1,-187.9,720. 

-37.1,-187.9,19. 

-37.1,-130.1,19. 

FGsDOOR #3 ON NEC. Y-SIDE 

4, F 

-37.1,-187.9,720. 

-5. ,-245.8,720. 

-5. ,-245.8,19. 

-37.1,-187.9,19. 

PG:WING SIMULATED PLATE #1 ON NEG. Y-SIDE 
3,F 

-109.35,-110.58,100. 

-120. ,-125. 192, 100. 

-120. ,-102.3928,-30. 

PG:W ING SIMULATED PLATE #4 ON NEG. Y-SIDE 

5, F 

-147. ,-106. 07,944.3 
-142.739,-468.34,867.57 
-135.083,-468.34,729.72 
-120. ,-188. ,458. 13 
-120. ,-118.84 ,458. 13 
PP: POLAR PLOT IN DB 
T 

1,2.766,3 

PD: ROLL PLANE CUT 

0 . ,0 . ,45 . 

0,360,1 
T, 50000. 

EX: EXECUTE 


The computed results are shown in Figure 108, however, measured 
results are not available. Comparing the two roll conical patterns, one 
can see that the blocking effect of the heat radiator causes large 
disturbances in the 30° to 60° and in the 100° to 160° range. This is 
due to the strong direct reflections from the heat radiator in the main 
beam region of the top mounted antenna. 
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(a) SIDE VIEW 



(cj TOP VIEW 


Figure 107. Computer simulated model for a crossed-slot antenna mounted 
on top of a Space Shuttle Orbiter when the payload bay 
doors are open. 
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Figure 108. Roll conical patterns (e p =45°) for a crossed-slot antenna 
mounted on top of a Space Shuttle Orbiter when the payload 
bay doors are open . 
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For the bottom mounted antenna case, a composite ellipsoid and 16 
plate model are used to simulate the basic structure. A computer 
simulated model and the resulting patterns are shown in Figures 109 and 
110, respectively. The input data is as follows: 


UN: IN INCHES 
3 

F6: SPACE SHUTTLE FOR THE BOTTOM MOUNTED ANTENNA WITH CLOSED PAYLOAD DOORS CASE 

40. . 108.. 2500.. 325. 

F 

0 . , 0 . , 0 . 

FC: 

T#F 

1055. . 0. 

FQ: FREQUENCY (2 GHZ) 

1 .2.. 1 . 

SG:CROSSED SLOT ( .39 *.78 WAVELENGTH) 

80., -4. 

2 

0 ., 0 . 

2.3.4.6.0. .1 .476.1 

1 .. 0 . 

0 ., 0 . 

2.3,4.6,90. ,1.476,1 
1 • 90 • 

PGjSIDE WALL SIMULATED PLATE #1 ON POS. Y-SIDE 

3, T 

-1 .,41 .5,-300. 

-155. ,109. ,-30. 

-5. ,107 .5,-30. 

PGjSIDE WALL SIMULATED PLATE #2 ON POS. Y-SIDE 

4, T 

-5.,1 07 .5,-30. 

-155. ,109. ,-30. 

-165. ,109.6, 18. 

-5. ,108. ,18. 

PG.-SIDE WALL SIMULATED PLATE #3 ON POS. Y-SIDE 
4»T 

-5. ,108. ,18. 

-165. ,109.6, 18. 

-165., 109.29,458. 

-5.,1 07 .69,458. 

PGjSIDE WALL SIMULATED PLATE #4 ON POS. Y-SIDE 

-5.,1 07 .69,458. 

-165., 109.29,458. 

-165. ,102.17, 1053. 

-5. ,100.57, 1033. 
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PG:WING SIMULATED PLATE #1 ON POS. Y-SIDE 

3, T 

1.. 106. 5.8. 

1.. 188.. 458. 

1.. 108.5.458. 

PG:WING SIMULATED PLATE #2 ON POS. Y-SIDE 
5,T 

1.. 108.5.458. 

1. . 188. .458. 

1.. 468.. 730. 

1.. 468. .867. 

1.. 107. .944. 

PG: BOTTOM SIMULATED PLATE #1 

4. T 

-1., -41 .5,-300. 

-155., -109., -30. 

-155. ,109. ,-30. 

-1., 41 .5,-300. 

PG: BOTTOM SIMULATED PLATE #2 
4.F 

-155., -109., -30. 

-165. .-109.6, 18. 

-165. ,109.6, 18. 

-155. ,109. ,-30. 

PG: BOTTOM SIMULATED PLATE #3 
4,F 

-165. ,-109.6, 18. 

-165. ,-109.29 ,458. 

-165., 109.29,458. 

-165. ,109.6, 18. 

PG: BOTTOM SIMULATED PLATE #4 
4,F 

-165. ,-109.29,458. 

-165. ,-102. 17, 1055. 

-165. ,102. 17, 1055. 

-165. ,109.29, 458. 

PG:SIDE WALL SIMULATED PLATE #1 ON NEG. Y-SIDE 

3, T 

-5., -107. 5, -30. 

-155. ,-109. ,-30. 

-1 .,-41 .5,-300. 

PG:SIDE WALL SIMULATED PLATE #2 ON NEG. Y-SIDE 

4, T 

-5. ,-108. ,18. 

-165. ,-109.6, 18. 

-155. ,-109. ,-30. 

-5., -107. 5, -30. 

PG:$IDE WALL SIMULATED PLATE #3 ON NEG. Y-SIDE 
4,T 

-5. ,-107 .69,458. 

-165. ,-109.29,458. 

-165. ,-109.6, 18. 

— 5. ,— 1 08. ,1 8. ?37 


PG:$IDE WALL SIMULATED PLATE #4 ON NEG. Y-SIDE 

4. T 

-5., -100.57, 1055. 

-165. ,-102. 17, 1055. 

-165. ,-109.29,458. 

-5., -107 .69,458. 

PG:WING SIMULATED PLATE #1 ON NEG. Y-SIDE 
3.T 

1., -108.5,458. 

1., -188., 458. 

1., -108 .5 ,8. 

PG:WING SIMULATED PLATE #2 ON NEG. Y-SIDE 

5. T 

1., -107., 944. 

1., -468., 867. 

1., -468., 730. 

1., -188., 458. 

1 .,-108.5,458. 

PP: POLAR PLOT IN DB 
T 

1 ,2.766,3 

PD: ROLL PLANE CUT 

0.,0.,45. 

0,360,1 
T, 50000. 

EX: EXECUTE 



(a) SIDE VIEW 



lb) FRONT VIEW 



(cj BOTTOM VIEW 


Figure 109. Computer simulated model for a crossed-slot antenna mounted 
on bottom of a Space Shuttle Orblter. 




<b) E # 


Figure 110. Roll conical patterns (9 p =45°l for a crossed-slot antenna 
mounted on bottom of a Space Shuttle Orbiter when the 
payload bay doors are closed . 
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When the payload doors are open, the whole structure Is still 
simulated by 16 plates as shown in Figure 111, although some of these 
plates are different from those obtained from the closed payload doors 
case. A roll conical pattern (e p =45°) is computed and then compared 
with that obtained with the payload doors are closed. The input data is 
as follows: 


UN: IN INCHES 
3 

FG: SPACE SHUTTLE FOR THE BOTTOM MOUNTED ANTENNA WITH OPEN PAYLOAD DOORS CASE 

40.. 108. .2500. .325. * 

F 

o.,o.,o. 

FC: 

T,F 

1055.. 0. 

FQ: FREQUENCY (2 GHZ) 

1 . 2 . . 1 . 

SG: CROSSED SL0T(.39 *.78 WAVELENGTH) 

80., -4. 

2 

0 . , 0 . 

2.3.4.6.0. . 1.476.1 

1 .. 0 . 

0 . ,0 . 

2.3.4.6.90. . 1.476.1 

1.. 90. 

FB: 

1 

6 

0.,41 .5,-300. 

0.,1 07 .5,-30. 

0., 107., 760. 

0., -107. ,760. 

0., -107. 5,-30. 

0.,-41 .5,-300. 

PGjSIDE WALL SIMULATED PLATE #1 ON POS. Y-SIDE 
3,T 

-1 .,41 .5,-300. 

-155. ,109. ,-30. 

-5., 107 .5,-30. 
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PGs SIDE MALL SIMULATED PLATE #2 ON POS. Y-SIDE 

5*T 

-5. ,107 .5,-30. 

-155. ,109. ,-30. 

-165. ,109 .29, 18. 

-80. ,108 .44,18. 

-5. ,107.69, 18. 

PGjSIDE MALL SIMULATED PLATE #3 ON POS. Y-SIDE 
4,T 

-5., 107 .69,1 8. 

-80. ,108 .44,18. 

-80., 1 08 .44,458. 

-5., 107 .69,458. 

PGsDOOR SIMULATED PLATE #1 ON POS. Y-SIDE 
4,F 

-80. ,108 .44, 18. 

—43.2,153.6,18. 

-43.2,153.6,458. 

-80., 108.44,458. 

PGsDOOR SIMULATED PLATE #2 ON POS. Y-SIDE 

4, F 

-43.2,153.6,18. 

-30.4,217 .6,18. 

-30.4,217.6,458. 

-43.2,153.6,458. 

PGsWING SIMULATED PLATE #1 ON POS. Y-SIDE 

3, T 

1.. 108.5.8. 

1.. 188. .458. 

1.. 108 .5.458. 

PGsWING SIMULATED PLATE #2 ON POS. Y-SIDE 

5, T 

1.. 108 .5.458. 

1.. 188. .458. 

1.. 468. .730. 

1 . ,468. ,867 . 

1.. 107. .944. 

PGs BOTTOM SIMULATED PLATE #1 

4, T 

-1.,-41 .5,-300. 

-155. ,-109. ,-30. 

-155. ,109. ,-30. 

-1 .,41 .5,-300. 

PGs BOTTOM SIMULATED PLATE #2 
4,F 

-155. ,-109. ,-30. 

-165. ,-109 .29, 18. 

-165. ,109 .29,18. 

109 .,-30. 
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PG: SIDE MALL SIMULATED PLATE #1 ON NEG. Y-SIDE 
3,T 

-5., -107. 5,-30. 

-155., -109., -30. 

-1 .,-41 .5,-300. 

PG:SIDE WALL SIMULATED PLATE #2 ON NEG. Y-SIDE 
5,T 

-5. ,-107.69, 18. 

-80.,-1 08.44,1 8. 

-165. ,-109.29,18. 

-155., -109., -30. 

-5^.-1 07.5,-30. 

PGjSIDE MALL SIMULATED PLATE #3 ON NEG. Y-SIDE 
4,T 

-5..-1 07 .69,458. 

-80., -108 .44 ,458. 

-80.,-1 08.44, 18. 

-5.,-107.69,18. 

PGjDOOR SIMULATED PLATE #1 ON NEG. Y-SIDE 

4.F 

-80., -108.44,458. 

-43.2,-153.6,458. 

-43.2,-153.6,18. 

-80. ,-108.44, 18. 

PGjDOOR SIMULATED PLATE §2 ON NEG. Y-SIDE 

4. F 

-43.2,-153.6,458. 

-30.4,-217.6,458. 

-30.4,-217.6,18. 

-43.2,-153.6,18. 

PGjMING SIMULATED PLATE #1 ON NEG. Y-SIDE 
3,T 

1., -108.5,458. 

1., -188. ,458. 

1., -108.5,8. 

PGjMING SIMULATED PLATE §2 ON NEG. Y-SIDE 

5, T 

1., -107., 944. 

1., -468., 867. 

1., -468., 730. 

1., -188., 458. 

1., -108.5,458. 

PPj POLAR PLOT IN DB 
T 

1,2.766,3 

PD j ROLL PLANE CUT 

0.,0.,45. 

0,360,1 

T,50000. 

EX j EXECUTE 
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The results shown in Figure 112 also indicate severe disturbances 
of the pattern due to the strong reflections from the payload doors. 
The above results indicate why one can lose use of his communication 
channel when the shuttle payload doors are open. 
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(a) SIDE VIEW 



(b) FRONT VIEW 



Figure 111. Computer simulated model for a crossed-slot antenna mounted 
on bottom of a Space Shuttle Orbiter when the payload bay 
doors are open. 
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Figure 112. Roll conical patterns (0p=45°) for a crossed-slot antenna 
mounted on bottom of a Space shuttle Orbiter when the 
payload bay doors are open . 
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